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CHAPTER  1 

INTRODUCTION 

The  study  of  flow  over  wavy  surfaces  has  attracted  attention  since  the 
work  of  Helmholtz  in  1868.  As  the  research  in  the  area  progressed  the  models 
became  more  sophisticated.  allowing  for  larger  amplitude  waves, 
incorporating  turbulent  stresses  and  interactions  between  the  turbulence  and 
the  wave. 

In  this  study  pressure  profiles  were  measured  over  a  wavy  surface  for 
turbulent  flow  in  a  rectangular  water  channel.  Using  a  voltage  response 
pressure  transducer  the  differential  pressure  was  measured  at  thirty-two 
locations  along  a  two  inch  wavelength.  All  pressures  were  referenced  to  the 
upstream  crest.  The  results  of  this  study  are  compared  to  several  other 
experimental  and  theoretical  works. 

The  motivation  for  this  research  is  threefold.  First,  it  allows  models  of 
turbulence  that  have  been  studied  over  simpler  geometries  (flow  over  a  flat 
plate,  for  instance)  to  be  tested  under  different  conditions.  This  gives  further 
insight  into  the  physics  of  turbulence.  The  second  reason  for  this  work  is  to 
gain  an  understanding  of  the  surface  forces  over  a  wavy  interface.  This 
information  can  be  used  in  studying  wave  generation.  The  third,  is  to  gain 
knowledge  of  the  effect  of  the  presence  of  the  wave  on  form  drag.  Previous 
studies  have  indicated  that  the  wave  can  reduce  skin  friction.  It  is  therefore 
necessary  to  know  the  increase  in  form  drag  to  determine  if  the  overall  drag  is 
actually  reduced  by  the  presence  of  the  wavy  surface. 
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CHAPTER  2 


LITERATURE  REVIEW 


The  first  model  for  flow  over  small  amplitude  waves  was  proposed  by 
Helmholtz  (1868).  This  neglected  viscous  effects  and  solved  the  equations  of 
motion  linearized  around  a  uniform  mean  velocity.  The  change  in  cross- 

sectional  area  of  flow,  due  to  the  presence  of  the  wave,  induces  a  variation  in 
the  mean  flow  along  the  wave  surface.  Due  to  the  Bernoulli  effect  this 
produces  a  pressure  variation  180°  out  of  phase  with  the  wave  surface. 

Miles  (1957)  dealt  with  the  generation  of  surface  waves  by  a  parallel 
shear  flow.  The  only  effect  of  turbulence  in  this  model  was  the  presence  of  a 
non-uniform  velocity  profile  u=u(y)  characteristic  of  a  turbulent  flow.  A 

mechanism  for  the  generation  of  waves  was  developed  using  the  inviscid  Orr- 

Sommerfeld  equation.  The  maximum  surface  pressure  was  calculated  to  shift 
from  the  trough  towards  the  downstream  crest.  Neither  the  Helmholtz  nor  the 
Miles  model  allowed  for  surface  shear  stresses  because  of  the  neglect  of 
viscosity. 

Benjamin  (1959)  and  Miles  (1959)  concurrently  reformulated  the 
problem  in  curvilinear  coordinates,  so  as  to  reduce  the  restrictions  on  the 
allowable  amplitude  of  the  perturbation.  These  models  also  include  the 

dominant  viscous  term  of  the  full  Orr-Sommerfeid  equation.  The  inclusion  of 
the  viscosity  was  found  to  have  an  effect  on  the  phase  relationship  of  the 
pressure  distribution.  These  models  do,  however,  neglect  wave-induced 

variations  of  turbulent  stresses  and,  therefore,  are  not  complete. 
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Thorsness  (1973)  also  worked  in  curvilinear  coordinates,  using  several 
turbulence  models  to  incorporate  wave-induced  turbulent  stresses.  Model  A  is 
basically  the  same  as  the  quasilaminar  approach  used  by  Benjamin.  Model  C 
accounts  for  Reynolds  stress  due  to  the  wave  through  a  turbulent  viscosity 
based  on  the  van  Driest  mixing  length  model  of  Loyd,  Moffat,  and  Kays  (1970). 
Model  D  also  used  mixing  length  considerations.  This  model,  however,  takes 
into  account  the  effect  of  the  acceleration  and  deceleration  of  the  flow  on  the 
turbulence  by  allowing  the  van  Driest  viscous  layer  thickness  parameter  to  be 
a  function  of  the  wall  pressure  gradient.  Model  D  was  found  to  give  better 
agreement  with  experiment  than  the  other  models. 

All  models  discussed  to  this  point  are  linear  and  limited  to  small 
amplitudes.  These  models  do  not  allow  for  separation  of  flow.  Abrams  (1979) 
used  a  nonlinear  boundary  layer  analysis  to  predict  occurrence  of  separation 
and  to  extend  the  range  of  calculations  to  larger  amplitude  waves.  He  also 
introduced  an  effective  pressure  gradient  into  Model  D  to  account  for  the  fact 
that  the  flow  cannot  adjust  to  the  pressure  gradient  instantaneously. 

In  further  work  Abrams  (1984)  observed  that  shear  stress  could  vary 
with  distance  from  the  wave  boundary  due  to  the  presence  of  pressure 
gradients.  For  this  reason  he  introduced  local  shear  stresses  into  Model  D, 
which  had  previously  recognized  only  wall  shear  stresses.  This  resulted  in 
Model  D*.  In  this  work  he  dealt  with  small  amplitude  waves  and  therefore  used 
linearized  equations  similar  to  Thorsness. 

Frederick  (1986)  developed  a  nonlinear  computer  code  for  use  with 
finite  amplitude  waves.  Using  spectral  methods  in  the  direction  of  flow  and 
finite  difference  methods  in  the  normal  direction,  the  code  solves  the 
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nonlinear  Reynolds-averaged  Navier-Stokes  equations.  The  code  modifies  the 
boundary  layer  program  of  McLean  (1983)  to  the  flow  geometry  of  waves  on 
the  lower  wall  of  a  rectangular  channel. 

Experimental  studies  of  pressure  variation  over  solid,  stationary, 
sinusoidal  surfaces  was  first  undertaken  by  Motzfeld  (1937)  in  a  wind  tunnel 
using  a  wavelength  of  30  cm  and  2a <j/k  =.05,  .10.  The  smaller  amplitude  wave 
gave  a  sinusoidal  pressure  variation  that  was  180°  out  of  phase  with  the  solid 
wave.  The  larger  amplitude  had  a  pressure  response  that  was  asymmetric,  the 
peak  shifted  slightly  downstream  from  the  trough. 

Zagustin,  Hsu,  Street,  and  Perry  (1966)  conducted  experiments  in  an 
open  water  channel.  Wall  surface  pressure  measured  over  2  and  3  foot 
wavelengths  with  2a<iA  =.04  and  .02,  respectively,  showed  linear  response. 

Kendall  (1970)  made  windtunnel  measurements  over  a  4  inch 
wavelength  with  2adM  =.0625.  The  wavy  surface  was  formed  by  a  smooth 

neoprene  rubber  sheet  which  was  mechanically  deformed  by  a  series  of  cams 
into  twelve  sinusoidal  waves.  These  waves  could  be  made  to  move  with  or 
against  the  airflow  at  controlled  speeds.  He  also  obtained  a  sinusoidal  pressure 
variation.  Saeger  and  Reynolds  (1971),  who  also  worked  with  moving  waves, 
obtained  a  sinusoidal  pressure  variation  for  2a<i/X  =.0094. 

Sigal  (1971)  made  measurements  in  a  wind  tunnel  over  aluminum 
sheeting  deformed  into  sinusoidal  waves.  Wavelengths  of  12  and  6  inches  with 
2a<j/X.  =.055  were  used,  with  five  waves  in  each  wavetrain.  The  pressure 
distributions  were  sinusoidal. 

Hsu  and  Kennedy  (1971)  made  measurements  of  airflow  in  a  wavy  walled 
pipe.  By  forming  fiberglas  over  wooden  mandrels  they  constructed  pipes 
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having  average  inner  diameter  of  .0408  feet  and  2ad/X  =.022  and  .011.  Direct 
comparison  of  results  is  questionable  here  due  to  the  differences  in  channel 
and  pipe  geometry. 

Lin,  Walsh,  and  Balasubramanian  (1984)  studied  drag  in  a  IS  inch  low 
turbulence  wind  tunnel.  Their  studies  included  a  variety  of  wave  forms.  Four 
of  these  were  sinusoidal  with  X=1  inch.  Values  of  2ad/X  were  0.01,  0.02,0.03,  and 
0.04.  A  0.S  inch  wavelength  with  2ad/X  =0.023  and  a  2  inch  wavelength  with 
2  ad  A  =0.005  were  also  used.  This  study  used  a  long  wavetrain  (35  waves)  so  as 
to  ensure  fully  developed  flow.  The  measurements  did  not  show  a  sinusoidal 
pressure  distribution  for  all  waves. 

In  this  lab  pressure  measurements  over  large  amplitude  waves  were 
done  by  Zilker  (1976),  Buckles  (1979),  and  Kuzan  (1983).  These  waves 
definitely  produced  large  separated  flow  and  are  therefore  not  comparable  to 
the  present  study. 

Cook  (1970)  studied  flow  over  waves  with  2a(j/X  =.05.  There  were  10 
waves  in  the  wavetrain.  This  wave  amplitude  is  likely  to  produce  separation  at 
lower  Reynolds  numbers.  Evidence  for  this  exists  in  the  nonlinearity  of  the 
pressure  profiles.  Studies  reported  by  the  other  investigators  with  2a<)/X  of 
this  value  or  higher  were  all  at  larger  Reynolds  numbers  and  thus 
nonseparated. 

Niederschulte  (1988)  conducted  experiments  under  almost  identical 
conditions  to  the  present  work.  In  the  current  study  measurements  were  made 
over  the  twenty-second  wave  of  a  twenty-four  wave  train.  Niederschulte's 
were  made  over  wave  eight  of  twelve.  For  2a j/X  =.03125  the  range  of 
dimensionless  wave  number,  a  +  =  2nv/Xu*,  was  the  same  for  both  studies. 
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Unfortunately  the  water  tunnel  employed  in  Niederschulte's  work  became 
unusable  due  to  age  and  deterioration  soon  after  measurements  over  a  wave 
with  2adA  =.0123  had  begun. 
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CHAPTER  3 


THEORY 


A.  Governing  Equations 

In  order  to  predict  turbulent  flow  properties  over  a  wavy  surface  of 
finite  amplitude  the  nonlinear  Reynolds-averaged  Navier-Stokes  equations  are 
considered.  In  cartesian  coordinates  the  two-dimensional,  steady  state 
momentum  equations  for  incompressible  flow  are  as  follows: 


1  aPd 
P  9x 


+  v 


32« 


32u 


3xz  3y2 


3x 


(  u'u') 


3y 


(-uV) 


1  3Pd 
P  3y 


+  V 


(3.1) 


(3.2) 


Here  p  is  the  density,  pd,  the  dimensional  surface  pressures,  the  kinematic 
viscosity,  u,  u'  and  v,  v',  the  time-averaged  and  fluctuating  velocities  in  the  x 
and  y  directions,  respectively.  The  turbulence  model  for  the  Reynolds  stress 
terms,  used  to  close  the  equations,  will  be  discussed  in  section  C. 

The  boundary  conditions  for  the  top  wall  and  the  wave  surface  are  no 
slip  and  no  penetration.  These  can  be  stated  as  follows: 

‘  u(x,  ad  cos  ax)  =  v(x,  a d  cos  ax)  =  0  ,  (3.3) 

«(x.yO  =  v(x,y,)  =  0,  (3.4) 

where  yt  |$  the  coordinate  of  the  top  wall  of  the  channel.  The  wave  surface  is 
described  by: 

y  =  ad  cos  ax,  (3.5) 
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with  aj  being  the  wave  amplitude,  a  =  2n/X  ,  the  wave  number,  and  X,  the  wave 
length.  In  the  horizontal  direction  the  periodicity  of  the  wave  surface  dictates 
the  following  boundary  conditions: 

u(0,y)  =  u(X„  y),  (3.6) 

v(0,y)  =  v(X„  y).  (3.7) 


B.  Numerical  Method  of  Frederick 

Linearized  theories  introduce  the  wavy  properties  of  the  flow  into  the 

« 

solution  by  proposing  that  flow  variables  consist  of  a  wave-averaged 
component  and  a  periodic  perturbation  term.  By  substituting  these  variables 
into  the  momentum  equations  and  eliminating  second  order  terms  in  the 
perturbation,  the  solution  of  problem  is  greatly  simplified.  This  imposes 
restrictions  on  the  amplitude  of  the  wave  that  are  most  likely  not  strictly  met 
in  the  present  study. 

Frederick  (1986)  dealt  with  this  problem  through  the  use  of  a  nonlinear 
channel  code.  The  code  is  used  to  solve  the  momentum  equations  (3.1)  and 
(3.2)  expressed  in  stream  function-vorticity  form.  The  stream  function,  y,  and 
the  vorticity,  o>,  are  defined  by 


u 


dy 

3y 


(3.8) 


v 


dx 


(3.9) 


3v  _  3u 
3x  dy 


(3.10) 
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In  order  Co  facilitate  computations  the  physical  domain  0  <=  x  <=  h  and 
a<j  cos  ax  <=  y  <=  yj  are  mapped  onto  the  rectangular  region  0  <=  e  <=  2n,  0  <=  q  <= 
Hi.  The  mapping  is  achieved  by  the  following  orthogonal  transformation: 


ax  =  e  + 


ay  =n  +  b, 


Zbi  f  cosh  Kq.-q)  ] 

_sin  ie  - — - 

.  «  L  s,nh  J 

•o  _  -i 

Zbi  sinh  i(nt-q) 

— cos  ie  - - 

l=L  ‘  L  ‘’ll  J 


(3-11) 

(3-12) 


The  coefficients  bj  are  chosen  to  approximate  the  wave  surface  at  q  =  0.  For  the 
waves  dealt  with  in  his  work  Frederick  found  eleven  terms  sufficient.  An 


additional  stretching  scheme  was  used  in  the  q  coordinate  to  give  better 
resolution .  near  the  upper  and  lower  boundaries,  where  the  velocity  gradients 
are  largest. 

The  vertical  derivatives  (those  with  respect  to  i)  )  were  evaluated  by 
standard  finite  difference  approximations.  The  horizontal  periodicity  of  the 
flow  field  allows  the  stream  function  and  vorticity  to  be  easily  represented  by 
a  Fourier  series.  Due  to  the  simplicity  of  calculating  Fourier  derivatives  the 
same  resolution  can  be  achieved  with  roughly  half  as  many  grid  points  as  are 
needed  with  a  finite  difference  scheme.  For  these  reasons  the  horizontal 


derivatives  were  evaluated  using  spectral  methods. 

Once  the  stream  function  was  calculated  the  cartesian  velocity  was 


recovered  by  the  chain  rule  equation: 

ft.  fr|!U»riS 

dy  oq  3y  de  3y 

For  more  detail  on  the  numerical  method  refer  to  Frederick  (1986). 


(3.13) 
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C.  Model  for  Turbulent  Stresses 

The  Reynolds  stress  terms  in  equations  (3.1)  and  (3.2)  can  be  modeled 
with  an  isotropic  eddy  viscosity: 


-uV  =  i-Rxx  =  vt2Sxx  =  vt2|H. 
P  dx 


-v  =TRyy=v«2syy=v«2lr 


3y 


-v'=F^-v'2s‘>-v-(!Hr) 


(3.14) 

(3.15) 

(3.16) 


where  Rjj  is  the  time-averaged  stress  tensor,  Sij  is  the  time-averaged  rate  of 
strain  tensor,  and  ut  is  the  turbulent  eddy  viscosity. 

Reynolds  and  Tiederman  (1967)  modeled  the  eddy  viscosity  as 


1  -exn 


fN72"5 


-y'r 


(3.17) 


where  h  is  the  average  half  channel  height,  K  is  the  von  Karman  constant,  y' 
is  distance  to  either  the  wave  surface  or  the  top  wall  (whichever  is  closer),  X  is 
the  local  shear  stress,  A  is  the  van  Driest  parameter  (a  measure  of  the 
thickness  of  the  viscous  wall  region). 

The  van  Driest  parameter  is  evaluated  using  Model  D*  developed  by 
Thorsness  (1975)  and  Abrams  (1984).  Model  D*  accounts  for  the  fact  that  the 
viscous  wall  region  thickens  with  increasing  negative  pressure  gradient  and 
thins  with  increasing  positive  pressure  gradient.  Loyd,  Moffat,  and  Kays 
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(1970)  suggest  the  following  functional  dependence  of  A  on  the  dimensionless 
pressure  gradient,  p+=  (dpd/dx)(v/pu*^): 

A  =  A  t  l  +  p*Kj  +p+  Kj  +  ...1  (3  18) 

where  A  is  the  flat  channel  value  of  A  and  Kj  and  K2  are  empirical  constants. 
For  the  nonequilibrium  condition  present  with  flow  over  waves  an  effective 
pressure  gradient,  p*eff  ,  should  be  used  in  (3.18).  The  effective  pressure 
gradient  is  defined  as  follows,  with  Kl  being  an  empirical  lag  constant: 

P+  —  P«vr 

(3l9> 

The  results  of  this  nonlinear  analysts  are  compared  with  the  experimental  data 
of  the  present  study. 
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CHAPTER  4 


EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


A.  Flow  System 

All  measurements  in  this  study  were  made  in  the  rectangular  water 
channel  designed  and  built  by  Niederschulte  (1988).  The  system  consists  of  an 
enclosed  horizontal  loop  approximately  50  feet  in  length  and  30  feet  in  height, 
through  which  water  circulates  in  a  counter-clockwise  direction,  as  pictured 
in  Figure  (4.1).  The  loop  has  a  capacity  of •1400  gallons  and  is  connected  to  two 
tanks,  one  of  1000  gallons  and  one  of  1650  gallons,  for  storage  of  test  fluids. 
The  system  was  constructed  using  316  or  304  stainless  steel,  plexiglas,  PVC, 
rubber,  or  rubber-lined  steel  to  avoid  corrosion. 

At  the  far  left  of  Figure  (4.1)  is  a  480  gallon  reservoir  tank  which 
contains  heating  and  cooling  coils,  providing  temperature  control  for  the 
channel.  With  careful  adjustment  of  the  heating/cooling  water  flow  rate  the 
temperature  of  the  loop  can  be  maintained  to  one-tenth  degree  centigrade.  A 
small,  one-half  horsepower  pump  can  be  used  to  circulate  the  contents  fo  this 
tank  through  a  bank  of  5pm  water  Filters. 

The  tank  empties  to  the  inlets  of  the  channel's  two  pumps,  which  are 
plumbed  in  parallel  and  are  run  only  one  at  a  time.  The  smaller  pump  is  a  five 
horsepower  Worthington  6CNG84  centrifugal,  316  stainless  steel  model  with  a 
six-inch  intake  and  six-inch  horizontal  discharge.  The  larger  pump  is  located 
in  a  lower  level  of  the  building  to  isolate  its  vibrations.  This  pump  is  a  316 
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Figure  4.1  Side  View  of  Flow  System 
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stainless  steel  Worthington  model  10S0-D  with  eight-inch  horizontal  intake 
and  six-inch  vertical  discharge.  It  has  a  60  horsepower  motor. 

The  parallel  pump  design  allows  for  a  wide  range  of  flowrates.  The 

small  pump  is  capable  of  operating  the  channel  through  the  range  of  40  g.p.m. 
to  800  g.p.m.  or  Reynolds  numbers  (based  on  bulk  velocity,  half-channel 

height,  and  water)  of  2400  to  43,000.  The  large  pump  can  deliver  2730  g.p.m. 
which  produces  a  Reynolds  number  of  160,000.  Both  pumps  discharge  to  tees 
in  the  channel/bypass  line.  The  bypass  line  leads  directly  back  to  the 
reservoir  tank  and  contaius  a  diaphragm  valve  to  control  its  flow  rate.  The 

channel  line  contains  a  butterfly  valve  to  control  flow  rate  to  the  channel. 
Manipulation  of  the  two  valves  allows  the  system  to  be  run  anywhere  from 
total  recirculation  (between  the  reservoir  tank  and  the  pump)  to  maximum 
channel  flow. 

The  channel  line,  made  of  eight-inch  PVC  piping,  contains  a  magnetic 

flowtube.  The  Foxboro  M/2800  series  flow  tube  measures  flow  using  the 
principles  of  magnetic  inductance.  The  advantages  of  this  form  of 

measurement  are  its  accuracy  (one-half  percent  of  full  scale  range)  and  low 
head  loss.  Another  major  advantage  is  the  fact  that  this  flow  measurement  is 
independent  of  temperature,  density,  and  viscosity  of  the  fluid  used.  The 
flowtube  outputs  its  reading  to  both  an  analog  meter  and  the  data  acquisition 
system  of  the  computer. 

The  flow  is  expanded  from  8  inch  diameter  to  22  x  30  inches  and  enters 
the  turning  elbows  of  the  system.  The  elbows  are  designed  to  decrease 
frictional  losses  and  undesirable  disturbances  resulting  from  abrupt  turning 
of  the  fluid.  A  honeycomb  is  used  to  straighten  the  flow  and  reduce  the 
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intensity  of  the  turbulence.  After  the  honeycomb  is  a  nozzle  system  that 

changes  the  cross  section  of  the  loop  from  22  x  30  inches  to  2  x  24  inches. 

The  channel  consists  of  four  70  inch  304  stainless  steel  segments,  each 
with  internal  cross  section  of  2  x  24  inches,  leading  to  the  test  section.  The  test 
section,  made  of  stainless  steel  and  plexiglas,  has  the  same  inner  dimensions  as 
the  channel  segments.  Deviations  from  stated  cross  sectional  dimensions  were 

limited  to  0.010  inches  for  the  first  three  channel  segments,  0.002  inches  for 
the  final  segment,  and  0.0002  inches  for  the  test  section.  (Test  section 
tolerance  excludes  plexiglas  components.)  The  bottom  wall  of  the  test  section 
consists  of  three  plexiglas  plates.  In  this  study  sinusoidal  waves  and  pressure 

taps  were  machined  onto  the  plates.  The  plates  will  be  discussed  in  greater 
detail  in  section  B. 

The  test  section  empties  into  the  reservoir  tank,  completing  the 

enclosed  flow  circuit.  The  design  criteria,  dimensions,  and  construction  of  the 
flow  system  are  described  in  detail  in  the  PhD  thesis  of  Nicderschulte  (1988). 

B.  Pressure  Apparatus 

The  equipment  used  to  measure  of  pressure  variation  over  a  wave 
consists  of  four  major  components:  plexiglas  plates  with  a  sinusoidal  surface 
and  pressure  taps,  a  pressure  transducer  and  demodulator,  a  header  device  to 
connect  the  taps  to  the  transducer,  and  a  device  to  calibrate  the  output  of  the 

transducer. 

The  plexiglas  plates  were  built  by  Bill  Knight  of  the  University  of 
Illinois  Chemical  Sciences  Machine  Shop.  The  plexiglas  slabs  were  first 
machined  flat  and  square.  Shoulders  for  mounting  the  plates  were  created  by 
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removing  material  from  the  perimeter  of  the  slab  to  leave  a  1  inch  wide 
border  of  1/2  inch  thickness.  (The  thickness  of  the  shoulder  was  increased  for 
strength  by  applying  3  inch  wide  by  1  inch  thick  plexiglas  pieces  to  the 
underside  of  the  finished  plate.)  Clearance  holes  for  1/4  inch  bolts  were 
drilled  1/2  inch  from  the  outer  edge  of  the  shoulder.  These  holes  were  spaced 

2  inches  apart.  Between  the  clearance  holes  and  the  body  of  the  slab  an  O-ring 

groove  was  milled  to  facilitate  sealing. 

The  sinusoidal  surfaces  were  formed  with  specially  milled  cutting  tools. 
Each  cutting  tool  was  formed  such  that  its  cutting  edge  was  a  complete  sine 
wave  with  an  additional  1/16  wave  period  on  each  end.  The  plates  were  set  up 

on  a  horizontal  milling  machine  and 'the  waves  were  cut  one  at  a  time  in  the 

spanwise  direction.  The  additional  length  of  the  cutting  tool  over  one  period 
allowed  the  machinist  to  overlap  the  waves,  reducing  machining  lines.  The 

0.03123  amplitude  waves  were  machined  onto  all  three  plates  of  the  test  section 
in  anticipation  of  use  in  other  research  projects.  The  0.0123  waves  were 
machined  only  onto  the  first  two  plates,  using  an  existing  flat  plate 

downstream.  The  downstream  plate  should  not  affect  the  pressure 

measurements,  which  are  done  over  the  second  plate. 

The  plates  were  polished  one  wave  at  a  time,  spanwise  across  the  entire 
length  of  the  wave.  Great  care  was  taken  to  remove  as  little  material  as 

possible  and  to  distribute  polishing  equally  to  crests  and  troughs.  Polishing 
was  accomplished  using  progressively  finer  agents:  300A,  600A  wet-dri 

sandpaper;  2/0,  3/0,  4/0  emery  paper;  and  finally  suspensions  of  aluminum 
oxide  in  deionized  water  (with  15p,  9.5p,  and  5p  particles).  The  final  finish  was 
glassiike  in  appearance. 
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The  pressure  taps  run  from  the  crest  of  the  twenty-second  wave  to  the 
crest  of  the  twenty-third  wave  in  either  a  twenty-four  wave  train  (amplitude  = 
0.0123  inches)  or  a  thirty-six  wave  train  (amplitude  =  0.031 25  inches).  The  taps 
are  0.040  inch  diameter  holes  drilled  through  the  wave  surface.  On  the  back 
surface  of  the  plate  the  taps  are  enlarged  to  0.127  inches.  This  enlargement 
allows  the  insertion  of  1/8  inch  stainless  steel  tubing  necessary  for 
connection  of  the  taps  to  the  pressure  transducer.  The  taps  meet  the  criteria 
outlined  by  Shaw  (1960)  for  minimizing  measurement  error. 

The  layout  of  the  pressure  taps  is  shown  in  Figure  (4.2).  There  are 
thirty-three  taps,  allowing  thirty-two  differential  pressure  measurements.  All 
measurements  were  referenced  to  the  upstream  crest.  The  taps  are  spaced  1/16 
inch  apart  in  the  streamwise  direction  and  1/8  inch  apart  in  the  spanwise 
direction.  The  arrangement  of  the  taps  (four  sets  of  eight,  each  set  at  an  angle 
of  43°  to  the  streamwise  direction)  allows  measurements  to  be  made  with  a 
streamwise  resolution  of  0.0623  inches,  yet  maintains  a  full  inch  distance 
between  taps  along  the  same  spanwise  location.  The  pressure  tap  layout  in 
centered  on  the  channel's  centerline. 

A  well  one  inch  deep  by  two  inches  wide  by  three  inches  long  was 
machined  into  the  back  of  the  waveplate,  centered  around  the  taps.  Stainless 
steel  tubing,  ranging  in  length  approximately  from  two  inches  to  ten  inches, 
was  inserted  into  each  tap  and  arranged  for  convenient  connection  to  the 
header  device.  The  well  was  filled  with  epoxy,  sealing  the  tubing  with  the 
plate  and  anchoring  the  tubing  arrangement.  A  short  length  of  tygon  tubing 
was  sealed  onto  each  piece  of  stainless  steel  tubing  with  General  Electric  RTV 
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Figure  4.2  Pressure  Tap  Layout 
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silicone  sealant.  The  tygon  tubing  is  used  to  connect  the  steel  tubing  to  the 
pressure  header. 

The  pressure  header  is  shown  in  Figure  (4.3).  This  was  built  by 

Niederschulte  and  modified  for  the  present  study.  It  was  built  with  1/8  inch 
stainless  steel  tubing  and  Swagelock  fittings  mounted  on  a  plexiglas  frame. 
The  header  consists  of  ten  lines  with  a  valve  in  each  line.  One  line  leads 

directly  to  the  low  pressure  inlet  of  the  transducer.  This  line  is  connected  to 
the  tap  at  the  upstream  crest  as  reference  for  all  measurements.  The  other 
nine  lines  are  connected  in  parallel  to  the  high  pressure  inlet  of  the 

transducer.  Eight  of  these  lines  are  used  to  connect  one  row  of  taps  to  the 
transducer.  The  ninth  line  can  be  used  as  a  bleed  line  to  remove  air  bubbles. 

Differential  pressure  was  measured  with  a  model  CJVR  pressure 

transducer  and  a  model  CJCD-4111  carrier  demodulator,  both  manufactured  by 
C.J.  Enterprises.  The  CJVR  has  a  full  scale  output  of  ±5  volts  and  depending  on 
the  pressure  diaphragm  used,  a  pressure  range  of  ±0.3  or±2.5  inches  of  water. 

The  pressure  transducer's  output  was  quite  sensitive  to  vibration  and 

position.  For  this  reason  the  transducer  was  anchored  to  the  header  device 
with  a  hose  clamp.  The  header  was  firmly  bolted  to  the  optical  traverse  used 

by  Niederschulte  (1988)  in  his  laser  doppler  experiments.  The  traverse  is  a 

large,  heavy  structure  designed  to  prevent  vibrational  effects  from  reaching 
that  which  is  mounted  on  it. 

A  calibration  device  was  constructed  (also  by  Niederschulte)  in  order  to 
obtain  the  imposed-pressure/voltage-output  relationship  for  the  pressure 
transducer  at  the  conditions  of  a  given  run.  The  calibrator  included  two  1/8 

inch  glass  tubes.  One  of  these  was  connected  via  the  header  to  the  reference 
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side  of  the  transducer,  the  other  was  connected  (also  via  the  header)  to  the 
measurement  side  of  the  transducer.  The  lines  were  filled  with  deionised 

water  and  air  bubbles  were  bled  out  through  bleed  screws  in  the  pressure 
transducer.  By  adding  water  to  one  tube  or  the  other  their  relative  height  of 
water,  and  thus  the  differential  pressure,  could  be  set  at  a  desired  level.  A  steel 
scale  graduated  in  1/64  inch  increments  was  mounted  next  to  each  tube.  The 
calibrator  was  leveled  with  a  bulls-eye  level  mounted  in  its  base.  The  body  of 
the  calibrator  was  precision  machined  stainless  steel  such  that  the  two  scales 
were  vertical  jthd  at  the  same  relative  height  if  the  bubble  in  the  bulls-eye 

level  w»j'  centered.  Therefore  the  imposed  head  could  be  accurately  measured. 
TK  calibration  device  is  shown  in  Figure  (4.4).  By  varying  the  head  and 

sampling  the  voltage  output  a  linear  calibration  relationship  was  found  for 

the  transducer. 

C.  Elflggdillfif 

Measurements  were  made  using  either  water  or  a  solution  of  dextrose  in 
water  as  the  operating  fluid.  The  use  of  fluids  of  different  viscosities  allowed 
measurements  to  be  made  over  a  wide  range  of  dimensionless  wave  numbers, 
a'*'  =  2xoAu*.  The  flow  loop  was  initially  filled  with  the  desired  fluid.  Before 
any  runs  were  made  with  that  fluid  the  flow  loop  was  run  at  various  speeds  for 
several  days  in  an  effort  to  remove  as  much  entrained  air  as  possible. 

Once  the  flow  rate  had  been  selected  for  a  given  run  the  loop  was 
configured  to  be  run  with  the  appropriate  pump.  Monitoring  the  flowrate 
with  the  magnetic  flowmeter,  the  bypass  and  butterfly  valves  were  adjusted 
until  the  desired  flow  was  achieved.  The  heating  or  cooling  water  was  run  to 
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bring  the  fluid  temperature  to  23°  C  and  then  adjusted  and  monitored  until  the 
fluid  held  that  temperature. 

Once  the  flow  system  was  set  the  pressure  taps  for  the  downstream  crest 
and  the  trough  were  connected  to  the  header  and  any  air  in  the  lines  was  bled. 
These  two  taps  are  where  the  extremes  in  the  pressure  profile  occur.  The  span 
adjust  and  zero  adjust  of  the  carrier  demodulator  were  adjusted  to  place  the 
pressure  range  occurring  over  the  wave  at  approximately  80%  of  full  scale. 

The  calibration  device  was  connected  to  the  header  and  any  air  bled  out. 
If  the  operating  lluid  was  dextrose,  all  dextrose  was  bled  out  also.  Fifteen 
different  heads  were  imposed  and  the  voltage  was  sampled. 

The  first  set  of  eight  pressure  taps  and  the  reference  tap  were 
connected  to  the  header,  and  bled  to  remove  air  (and  water  introduced  by  the 
calibration,  if  the  operating  fluid  was  dextrose).  The  pressure  at  each  tap 
location  was  measured  by  opening  the  valve  on  its  connecting  line  and 
sampling  the  pressure  transducer's  voltage  output.  A  twelve  bit  analog  to 
digital  converter  manufactured  by  Metra-Byte  Corporation  was  used  to  collect 
the  data  on  an  IBM  AT  Personal  Computer.  Spectral  analysis  of  several  data  sets 
indicated  that  random  noise  overpowered  signal  at  approximately  30  Hz.  Stable 
averages  were  obtained  with  sample  times  of  100  seconds.  The  voltages  were 
sampled  at  100  Hz  for  130  seconds.  The  voltages  were  also  monitored  with  a 
digital  voltage  meter.  In  this  manner  erroneous  voltages  due  to  air  bubbles  in 
the  lines  could  be  detected,  the  line  bled,  and  the  point  retaken.  After  all  eight 
taps  were  sampled  the  next  set  was  connected  and  the  process  repeated.  After 
all  four  sets  were  finished  a  second  calibration  was  conducted. 
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The  relationship  between  imposed  differential  pressure  and  voltage 
output  was  linear.  The  correlation  coefficient  was  always  above  0.980.  The 
deviation  from  a  correlation  of  1.00  was  due  to  scatter  in  the  points  as  opposed 
to  any  curvature  in  the  relationship  of  imposed  differential  pressure  versus 
voltage  output.  The  two  calibrations  agreed  within  two  percent  of  one  another 
with  regard  to  the  slope  of  the  plot.  The  intercept  of  the  line  often  changed 
from  one  calibration  to  the  other.  It  was  observed  that  the  excessive  bleeding 
necessary  for  the  connection  of  the  calibration  device  sometimes  caused  a 
shift  in  the  output  voltage.  These  shifts  were  not  noticed  with  the  changing  of 
the  pressure  taps  within  a  run  as  they  require  considerably  less  bleeding  of 
the  lines.  As  a  result  the  intercept  chosen  for  analyzing  the  data  was  the  one 
that  would  position  the  pressure  profile  so  that  the  differential  pressures 
seemed  most  reasonable.  It  is  important  to  note  that  the  intercept  only  affects 
the  zeroth  order  harmonic,  or  DC  component.  This  component  has  no  effect  on 
the  quantities  calculated  in  this  study,  such  as  phase  angles  and  amplitudes  of 
the  profile,  or  drag  coefficients.  A  sample  calibration  is  shown  in  Appendix  C. 

Data  was  stored  in  binary  form  on  the  hard  disc  of  the  IBM.  A  Fortran 
program  (Press.For,  see  Appendix  B)  on  the  IBM  was  used  to  unpack  the  data  to 
standard  ASCII  form  as  voltages.  The  program  then  converted  the  voltages  to 
pressure  with  the  linear  relation  obtained  from  calibration.  The  pressures 
were  non-dimensionalized  by  dividing  through  by  pu*2,  where  u*  is  the 
friction  velocity  (equal  to  the  square  root  of  the  wall  shear  stress  divided  by 
density).  The  mean  pressure  gradient  (that  which  would  be  present  in  a  flat 
channel)  was  removed  by  assuming  a  linear  pressure  drop  across  the  two  inch 
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wavelength.  The  magnitude  of  the  drop  was  assumed  to  be  described  by  the 
relation 


v 

h  u  * 


(4.1) 


from  the  correlations  of  Dean(1978).  Plotting  files  of  non-dimensionalized 
pressure  versus  location  on  the  wave  were  the  output  of  this  program.  The 
plotting  files  were  networked  to  a  Digital  Equipment  Company  microVax.  On 
the  Vax  another  Fortran  program  (Fit3u.for,  see  Appendix  B)  fit  the  pressure 
profiles  to  three  harmonics  through  Fourier  analysis.  The  harmonic  files  and 
the  data  files  were  accessed  by  an  Apple  Macintosh  computer  for  plotting. 
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CHAPTERS 

RESULTS 

Differential  pressure  profiles  were  measured  for  eight  different 
Reynolds  numbers  over  the  0.03 125  inch  amplitude  wave  and  for  ten  Reynolds 
numbers  over  the  0.0123  inch  amplitude  wave.  Figures  (5. 1-5.8)  show  the 
profiles  measured  over  the  0.03125  inch  wave  and  figures  (5.9-5.18)  show  the 
profiles  measured  over  the  0.0123  inch  wave.  These  figures  show  differential 
pressure  (p  -  Pcrest)  nondimensionalized  with  fluid  density,  p,  and  the  square 
of  the  friction  velocity,  u*.  plotted  versus  streamwise  location  along  the  wave, 
nondimensionalized  with  the  wavelength,  X.  The  friction  velocity  is  defined  as 

u*=ub  (Cf/2)0-5  (5.1) 

where  u(,  is  the  bulk  velocity  and  Cf  is  the  friction  factor.  The  friction  factor 
is  obtained  through  the  correlation  of  Dean(1978)  for  flow  through  a 
rectangular  channel  (this  correlation  is  based  on  smooth  walls) 

Cf  =  0.0614  Re’0  25  (5.2) 

where  Re  is  the  Reynolds  number  based  on  half  channel  height.  The 
experimental  data  and  run  conditions  are  tabulated  in  Appendix  A. 

The  pressure  profiles  were  found  to  depend  on  the  dimensionless  wave 
number,  a+=  2xv/X.u*,  and  the  dimensionless  amplitude  of  the  wave  surface,  a+= 
a(ju*/v.  The  pressure  response  can  be  characterized  as 

p(x/A.)  =  2]  a;  cos(2xx/X  +  9j)  (5.3) 

where  9i  is  the  phase  angle  between  the  pressure  response  and  the  wave 
surface.  A  positive  phase  angle  indicates  the  pressure  maximum  is  shifted 
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upstream  from  the  wave  crest.  The  profiles  taken  at  the  three  smallest  a+  for 
the  0.03 125  inch  wave  have  a  linear  response  and  are  adequately  described  by 
one  term  in  equation  (3.1).  The  rest  are  adequately  described  by  three  terms. 
Over  the  0.0123  inch  wave  all  the  profiles  except  that  for  the  highest  a+  are 
adequately  described  by  one  term.  In  order  to  compare  the  ai  of  runs  using 
working  fluids  of  different  kinematic  viscosity  it  is  necessary  to  look  at  the 
ratio  ai/a+.  (The  relationship  of  ai  and  0j  to  o+  and  a+  will  be  covered  in 
Chapter  6.) 

Figures  (3. 1-3.4)  show  the  profiles  taken  over  the  0.03123  inch  wave 
using  water  as  the  working  fluid.  Figure  (3.1)  was  taken  at  a  Reynolds 
number  of  41,600  and  a+  of  0.0016.  The  ratio  ai/a+  is  0.33  and  the  phase  0}  is 
167®.  This  profile  is  described  quite  well  by  a  linear  response,  the  amplitude  of 
the  second  harmonic  being  three  percent  of  the  first.  This  value  is  due  more 
to  experimental  scatter  the  any  distortion  of  the  profile.  The  curve  created 
using  three  harmonics  falls  directly  over  that  using  a  single  harmonic  term. 

Figure  (3.2)  is  the  profile  taken  at  a  Reynolds  number  of  30,400  and  a+ 
of  0.0021.  The  ratio  ai/a+  is  0.625  and  9i  is  171°.  This  profile  is  also  linear  with 
*2/*l  equal  to  0.03.  The  three  harmonic  curve  falls  almost  exactly  over  the 
one  harmonic  curve. 

Figure  (5.3)  was  obtained  at  a  Reynolds  number  of  21,600  and  a+  of 
0.0029.  The  ratio  ai/a+  is  0.77  and  9i  is  168°.The  second  harmonic  amplitude  is 
six  percent  of  the  first.  In  this  case  the  three  harmonic  curve  diverges  from 
the  single  harmonic  curve,  particularly  near  the  crests  of  the  wave  surface. 
This  is  considered  a  borderline  linear/nonlinear  profile.  Note  that,  in  this 
profile  (as  well  as  the  other  profiles  to  come),  although  the  additional 
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harmonic  terms  may  take  on  significant  amplitudes  their  relative  phases  are 
such  that  they  affect  the  shape  of  the  profile  withont  causing  significant 
change  to  the  magnitude  of  the  profile.  In  other  words,  even  though  the 
curve  created  using  three  harmonic  terms  differs  from  the  single  harmonic 
curve,  the  first  amplitude  term  a|,  adequately  describes  the  magnitude  of  the 
pressure  response. 

Figure  (5.4)  is  the  profile  taken  at  a  Reynolds  number  of  14,900  and  a+ 
of  0.0040.  The  amplitude  ratio  aj/a+  is  0.82  and  0j  is  169°.  This  profile  is 
nonlinear  with  a2/ai  equal  to  0.11.  Although  nonlinearity  has  set  in  the 
phase  is  still  similar  to  that  of  the  linear  profiles. 

Figures  (5. 5-5. 8)  show  the  profiles  taken  over  the  0.03125  inch  wave 
using  dextrose  solution  as  the  working  fluid.  Figure  (5.5)  is  the  profile  taken 
at  a  Reynolds  number  4690  and  a+  of  0.0108.  The  amplitude  ratio  ai/a+  is  1.40 
and  6j  is  150°.  The  profile  is  quite  nonlinear  with  a2/aj  equal  to  0.14.  The 
phase  shift  is  now  significantly  lower  than  that  observed  for  the  smaller  a+ 
runs. 

Figure  (5.6)  shows  the  profile  measured  at  a  Reynolds  number  of  3250 
and  a+  equal  to  0.0150.  The  ratio  aj  /a+  is  1.49  and  9  j  is  149°.  The  second 
harmonic  amplitude  is  eighteen  percent  of  the  first  harmonic's.  The  figure 
(5.7)  is  the  profile  taken  at  a  Reynolds  number  of  3220  and  a+  equal  to  0.0151. 
This  run  was  made  to  establish  the  level  of  reproducibility  of  the  experiments. 
The  ratio  ai /a+  is  1.57  and  0}  is  144°.  The  second  harmonic's  amplitude  is 
fourteen  percent  of  the  first's.  The  ratio  ai  /a+  agrees  to  5.1%  and  the  phase 
angle  agrees  to  5°.  Both  profiles  are  nonlinear  and  have  similar  shapes. 
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Figure  (5.8)  shows  the  profile  taken  at  a  Reynolds  number  of  2360  and 
a*  equal  to  0.0198.  The  ratio  ai/a+  is  1.36  and  6i  is  146*.  The  ratio  a2/ai  is  0.20 
and  the  profile  is  highly  nonlinear. 

Figures  (3.9-3.12)  show  the  profiles  obtained  over  the  0.0123  inch  wave 
using  water  as  the  working  fluid.  Figure  (3.9)  is  the  profile  from  the  run 

taken  at  a  Reynolds  number  of  66,400  and  a4*  equal  to  0.0011.  The  ratio  ai/a+  is 

0.432  and  0|  is  173°.  The  ratio  a2/ai  is  0.06.  This  is  probably  inflated  due  to 

experimental  scatter  as  the  three  harmonic  profile  follows  the  one  harmonic 

profile  very  closely. 

Figure  (3.10)  was  measured  at  a  Reynolds  number  of  43,900  and  a+  equal 
to  0.00134.  The  amplitude  ratio  aj/a+  is  0.302  and  0|  is  172°.  The  amplitude  of 
the  second  harmonic  is  four  percent  of  the  first  and  the  profile  is  very  linear. 

Figure  (3.11)  is  the  profile  measured  at  a  Reynolds  number  of  32,100 
and  a+  equal  to  0.0020.  The  ratio  ai/a+  is  0.368  and  0j  is  172°.  The  amplitude  of 
the  second  harmonic  is  again  four  percent  of  the  first  and  the  profile  is 
linear. 

Figure  (3.12)  shows  the  pressure  profile  taken  at  a  Reynolds  number  of 
21,400  and  a+  equal  to  0.0029.  The  value  of  ai/a+  is  0.670  and  ©i  is  172°.  The 

ratio  a2/ai  is  0.08,  but  this  is  due  mostly  to  experimental  scatter  as  the  profile 

is  linear. 

Figures  (3.13-3.18)  show  the  profiles  taken  over  the  0.0123  inch  wave 
using  dextrose  solution  as  the  working  fluid.  Figure  (3.13)  shows  the  pressure 

profile  obtained  at  a  Reynolds  number  of  13,000  and  a+  equal  to  0.0039.  The 

ratio  ai/a+  is  0.993  and  0  (  is  163°.  The  second  harmonic  amplitude  is  six 
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percent  that  of  the  first  and  the  profile  is  linear.  The  profile  phase  has  begun 
to  shift  from  the  low  a+  values. 

Figure  (5.14)  was  taken  at  a  Reynolds  number  of  9450  and  a*  equal  to 
0.0059.  The  ratio  ai/a+  is  1.033  and  0i  is  162°.  The  value  of  a2/aj  is  0.05  and  the 
pressure  response  is  linear. 

Figure  (5.15)  the  profile  measured  at  a  Reynolds  number  of  6730  and  a+ 
equal  to  0.0079.  The  ratio  aj/a +  is  1.172  and  ©i  is  157°.  The  ratio  of  a2/aj  is  0.09 
and  the  profile  is  a  borderline  linear/nonlinear  pressure  response. 

Figure  (5.16)  shows  the  measured  profile  at  Reynolds  number  5200  and 

a+  equal  to  0.0099.  The  ratio  aj/a+  is  1.203  and  0 ]  is  162°.  This  is  another 

borderline  profile  with  the  amplitude  of  the  second  harmonic  being  eight 
percent  of  the  first.  Figure  (5.17)  is  the  profile  taken  at  a  Reynolds  number  of 
5120  and  a+  equal  to  0.0101.  The  ratio  a|/a+  is  1.159  and  01  is  156°.  The  second 
harmonic  has  an  amplitude  that  is  six  percent  that  of  the  first  and  the  profile 
is  also  borderline.  The  amplitudes,  aj ,  agree  to  three  percent.  The  variation  is 
0]  is  perhaps  due  to  one  or  two  erroneously  low  points  over  the  0.4  to  0.5  xA 
range  in  figure  (5.16). 

Figure  (5.18)  shows  the  pressure  profile  measured  at  a  Reynolds  number 
of  3300  and  a+  equal  to  0.0148.  The  ratio  ai/a+  is  1.420  and  0t  is  164°.  This  is  a 
larger  than  expected  phase  angle.  One  explanation  is  that  with  a2/ai  equal  to 
0.12  this  is  the  first  truly  nonlinear  profile  measured  over  the  0.123  inch 
wave.  The  three  harmonic  curve  (which  follows  the  data  well),  has  its  maxima 

located  to  the  right  of  the  one  harmonic  maxima.  This  indicates  the  true  phase 

angle  is  lower  than  0i . 
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It  should  be  noted  that  some  of  the  profiles  contain  an  erroneous  zeroth 
order  harmonic,  shifting  the  profile  up  or  down  by  a  constant  amount.  This  is 
due  to  the  calibration  difficulties  discussed  in  the  chapter  on  Experimental 
Equipment  and  Procedures.  This  shift  does  not  affect  the  amplitudes  or  phases 
of  the  profile  nor  does  it  affect  the  form  drag  coefficients  calculated  in  the 
next  chapter. 
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Figure  5.1  Pressure  Profile  2a <j/X  =  0.03125,  a+=  0.0016 
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Figure  5.2  Pressure  Profile  2adA  =  0.03125,  a+=  0.0021 


Pressure  Profile  2nd A,  =  0.03125,  a+=  0.0029 


Figure  5.4  Pressure  Profile  2ndA  =  0.03125,  a+=  0.0040 
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Figure  5.6  Pressure  Profile  2adA  =  0.03125,  a+=  0.015 


Figure  5.8  Pressure  Profile  2a<jA  =  0.03125,  a+=  0.020 
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Figure  5.12  Pressure  Profile  2a<jA  =  0.0123,  a+=  0.0029 
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Figure  5.13  Pressure  Profile  2adA  =  0.0123,  a+=  0.0039 
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Figure  5.14  Pressure  Profile  2adA  =  0.0123,  a+=  0.0059 
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Figure  5.15  Pressure  Profile  2a d/X  ~  0.0123,  a+=  0.0079 
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Figure  5.16  Pressure  Profile  2a<jA  =  0.0123,  a+=  0.0099 
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Figure  5.17  Pressure  Profile  2adA  =  0.0123,  a+=  0.0101 
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Figure  5.18  Pressure  Profile  2adA  =  0.0123,  a+=  0.0148 
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CHAPTER  6 

DISCUSSION  OF  RESULTS 

The  measurements  made  in  this  study  covered  a  wider  range  of  the 

dimensionless  wave  number,  a + ,  than  any  of  the  previous  experimental 
studies.  Because  of  this  they  are  of  great  use  in  evaluating  the  various 

theoretical  works  discussed  in  the  literature  survey.  The  pressure  response  to 
flow  over  a  wavy  surface  can  be  summarized  by  the  ratio  of  the  first  harmonic 
amplitude  to  the  dimensionless  wave  amplitude,  ai/a+,  and  the  phase  angle  of 
the  first  harmonic,  9j .  At  larger  dimensionless  wave  numbers  the  shape  of  the 
pressure  response  is  markedly  distorted  from  the  simple  one  harmonic 

approximation.  By  comparing  profiles  taken  over  the  same  wave  at  a+  that  are 
almost  the  same,  the  level  of  reproducibility  can  be  seen  as  roughly  five 

percent  for  first  harmonic  amplitude,  ai,  and  five  degrees  for  the  phase  angle, 
0 1 .  The  variation  in  magnitude  and  phase  of  the  actual  measured  profile  from 
the  linear,  one  harmonic  approximation  is  within  these  tolerances  (for  all  but 

the  highest  o+)  .  This  is  not  to  say  that  these  profiles  are  exactly  linear,  but 

only  that  aj/a+  and  6i  are  appropriate  for  comparison. 

The  data  from  the  work  of  Lin  et  al  (1984),  Niederschulte(1988),  and  this 
study  are  summarized  in  Table  (6.1).  Figure  (6.1)  compares  the  amplitude  ratio, 
*l/s+,  of  these  experimental  studies  with  linear  predictions  from 
Thorsness(1975)  using  Model  A,  Abrams(1984)  using  Model  D*,  and  the 

nonlinear  analysis  of  Frederick!  1986).  Figure  (6.2)  compares  phase  angles,  01, 
in  a  similar  fashion.  In  both  of  the  figures  and  the  table,  the  work  of  Lin  et  al 
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Table  6.1  SUMMARY  OF  EXPERIMENTAL  DATA 
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Table  6.1  SUMMARY  OF  EXPERIMENTAL  DATA,  continued 
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Table  6.1  SUMMARY  OF  EXPERIMENTAL  DATA,  continued 


alpha  + 

Figure  6.2  Comparison  of  Phase  Angles 
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is  referenced  as  NASA  and  the  work  of  Niederschulte  is  referenced  as  Nieder. 

Lin  et  al  measured  pressure  over  waves  with  similar  2ad/X  to  the 
present  work.  They  cover  an  a+  range  of  0.0027  to  0.0063.  The  amplitude  ratios 
are  in  good  agreement  with  the  present  study,  generally  within  the 
reproducibility  of  the  work.  The  phase  angles  agree  for  a+  above  0.0022. 
Below  that,  however,  the  NASA  data  indicates  a  faster  approach  to  the  expected 
asymptote  of  0]=  180°.  A  possible  cause  for  this  disagreement  is  the  different 
dimensions  of  the  flow  system  used  in  the  two  studies.  The  NASA  data  was 
obtained  in  a  wind  tunnel  where  the  top  and  bottom  walls  were  fifteen  inches 
apart.  The  water  channel  in  this  study  has  a  distance  of  two  inches  between 
the  top  and  bottom  walls,  which  can  be  compared  to  a  boundary  layer 
thickness. 

The  data  of  Niederschulte  falls  into  two  regions:  The  first  was  obtained 
with  water  over  a+  of  0.00191  to  0.0022;  the  second  was  obtained  with  dextrose 
solution  over  a+  of  0.00913  to  0.02245.  The  data  taken  in  the  lower  a+  region 
agree  well  with  the  present  results  in  amplitude  ratio,  ai/a+.  The  phase  angles 
are  somewhat  lower  (approximately  four  percent)  than  those  obtained  in  this 
study.  The  phase  angle  measured  at  a+  equal  to  0.0016  over  the  0.03125  inch 
wave  in  this  study  does  agree  with  Niederschulte's  results.though. 

The  phase  angles  from  Niederschulte's  high  a*  measurements  agree 
quite  well  with  those  of  the  present  work.  A  comparison  of  the  shape  of 
profiles  from  the  two  studies  at  similar  a  +  shows  remarkable  similarity 
considering  the  data  were  collected  in  two  different  flow  systems,  using 
different  wave  plates  and  differing  arrangements  for  connecting  the  pressure 
taps  to  the  pressure  transducer. 
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The  amplitudes  of  Niederschulte's  high  a+  data  present  some  difficulty. 
The  values  of  ai/a+  are  a  much  as  fifteen  percent  higher  as  measured  by 
Niederschulte.  The  amplitude  ratio  also  reaches  a  maxima  at  a+  equal  to  0.0171. 
The  amplitude  ratio  measured  over  the  0.0123  inch  wave  does  not  agree  with 
the  0.03125  inch  wave  data. 

It  appears  that  Niederschulte's  amplitude  data  are  suspect.  The  first 
argument  supporting  this  statement  is  the  abruptness  of  the  drop  in  amplitude 
after  the  maxima  is  reached.  Also,  the  values  at  the  two  highest  a+  agree  with 
the  current  study,  although  that  may  only  be  coincidence.  Secondly,  the  data 
point  taken  over  the  0.0123  inch  wave  does  not  agree  with  his  other  data,  while 
the  data  from  the  two  wave  amplitudes  agree  quite  well  in  the  present  study. 
Finally  the  data  of  the  current  study  follow  a  recognizable  trend  over  a  wide 
range  of  a+.  The  low  a+  data  of  Niederschulte  are  effectively  clumped  into  two 
data  points,  making  trends  very  difficult  to  discern.  The  high  a+  data  shows 
very  little  pattern.  Based  on  the  similarity  of  Niederschulte's  profile  shapes  to 
the  current  study,  it  may  be  concluded  that  the  profiles  were  measured 
accurately,  but  the  pressure  transducer's  voltage  output  to  imposed 
differential  pressure  relationship  was  incorrectly  calibrated  throughout  the 

high  a*  runs.  These  runs  were  all  made  using  dextrose  solution,  which  makes 
calibration  more  difficult  as  the  need  to  bleed  the  lines  of  the  pressure  header 
is  increased  due  to  the  presence  of  the  dextrose. 

The  phase  angles  of  this  study  are  best  represented,  among  the  models 
considered,  by  Abrams'  linear  calculations  using  Model  D*.  All  of  the 

experimental  data  presented  here  follow  the  shape  of  Abrams’  curve  but  are 

shifted  a  few  degrees  lower  than  his  predictions.  The  notable  exceptions  are 
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the  phase  angles  at  a+  equal  to  0.0099  and  0.0148  over  the  0.0123  inch  wave, 
which  are  significantly  higher.  The  phase  angle  at  a+  equal  to  0.0101  (which 
should  be  the  same  as  that  of  a+  equal  to  0.0099)  falls  where  expected  from  the 
trend  of  the  data. 

The  amplitude  ratio  predictions  of  Frederick's  channel  simulation  are  in 
excellent  agreement  with  the  experimental  data  of  this  study.  The  NASA  data 
appear  in  a  range  where  all  three  theoretical  models  predict  similar  results. 

The  data  of  Niederschulte  don't  agree  with  any  of  the  models,  but  this  has 
already  been  discussed. 

As  stated  in  the  introduction,  one  motivation  for  researching  flow  over 
a  wavy  surface  is  the  reduction  of  drag.  Skin  friction  can  be  reduced  by 
imposing  a  wavy  boundary  on  the  flow.  The  wavy  surface,  however,  will 

inevitably  increase  form  drag.  If  the  increase  in  form  drag  is  larger  than  the 

decrease  in  skin  friction  then  no  gain  is  made. 

It  is  of  interest,  therefore,  to  see  if  a  drag  coefficient  calculated  using 
the  one  harmonic  approximation  to  the  pressure  profile  gives  an  accurate 
representation  of  the  form  drag  or  if  it  is  necessary  to  use  the  actual  pressure 
profile  in  the  calculation.  A  nondimensional  drag  coefficient  was  calculated 
by  dividing  the  form  drag  the  wave  exerts  on  the  fluid  by  the  effective  area 
(the  area  present  in  the  absence  of  a  wave,  equal  to  the  wavelength  times  the 
width)  and  pu*2. 

Cp  =  lApu*2.J  pd  sin  0  dl  (6.1) 

where  f)  =  arctan(dy/dx  of  the  wave  surface)  and  dl  is  the  differential  distance 

along  the  wave  surface.  The  drag  coefficients  were  calculated  using  both  the 
one  harmonic  fit  of  the  pressure  profile  (a  linear  approximation  to  the 
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measured  profile)  and  the  three  harmonic  fit  (which  follows  the  data  very 
closely). 

The  calculated  drag  coefficients  are  plotted  in  Figure  (6.3).  The 
coefficients  for  the  0.0123  inch  wave  show  close  agreement  between  those 
calculated  using  the  one  harmonic  approximation  and  those  using  the  more 
accurate  three  harmonic  fit.  The  agreement  is  not  as  close  for  the  0.03123 
inch  wave  were  the  two  calculations  generally  differ  by  fifteen  percent.  This 
is  reasonable  since  the  one  harmonic  fit  and  the  three  harmonic  fit  tend  to  be 
quite  similar  for  the  0.0123  wave  and  not  as  similar  for  the  0.03123  wave.  The 
one  harmonic  calculation  may  still  be  appropriate  for  the  0.03123  inch  wave  as 
the  more  exact  calculation  gives  a  larger  coefficient  for  some  of  the  cases  and 
a  smaller  coefficient  in  the  rest  of  the  cases.  Also  included  in  the  graph  is  a 
form  drag  coefficient  calculated  by  Frederick.  It  is  in  good  agreement  with 
those  of  this  study. 

Another  important  consideration  in  flow  over  a  wavy  surface  is  the 
possibility  of  flow  separation,  or  areas  of  stagnant  flow  on  the  downward 
sloping  face  of  the  wave.  Separated  flow  is  predicted  when  the  surface  shear 
stress  is  allowed  to  reach  zero.  In  his  nonlinear  analysis  Abrams  (1979) 
mapped  the  boundary  between  separated  and  nonseparated  flow  as  2a<jA 
versus  <x+.  Figure  (6.4)  shows  where  the  data  of  this  study  lie  on  that  map.  As 
seen  on  the  map  all  of  the  flows  studied  over  the  0.0123  inch  amplitude  wave 
should  be  nonseparated.  The  pressure  profiles  support  this  as  none  of  them 
are  truly  distorted.  The  map  predicts  that  flow  over  the  0.03125  inch  amplitude 
wave  should  separate  above  an  a+  of  approximately  0.006.  The  data  support 
this  prediction  also.  The  profiles  for  a+  below  0.003  are  quite  symmetric.  The 


Figure  6.3  Calculated  Drag  Coefficients  versus  Alpha 
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profile  at  a+  equal  to  0.004  begins  to  show  some  distortion.  The  rest  of  the 
profiles  are  at  a+  greater  than  0.01  and  are  all  distorted.  This  distortion  could 
be  due  to  a  separated  zone. 

In  conclusion,  it  is  in  the  high  a'*'  regions  that  linear  approximations 
should  break  down.  In  this  region  the  nonlinear  predictions  of  amplitude  by 
Frederick  match  the  experimental  data  best.  Surprisingly,  the  predictions  of 
phase  angle  that  match  the  experimental  data  best  are  from  the  linear 
calculations  using  Model  D*  by  Abrams. 

Although  it  appears  that  the  amplitude  of  the  surface  waves  in  this 

study  are  large  enough  to  necessitate  nonlinear  calculations,  the  pressure 
response  seems  to  be  satisfactorily  described  as  a  first  approximation  with  a 

linear  (one  harmonic)  approach.  This  statement  can  be  supported  by  the 

observance  that  the  pressure  response's  amplitude  and  phase  as  reported  by 

the  one  harmonic  approximation  are  within  the  experimental  reproducibility 
of  the  measurements  to  the  actual  amplitude  and  phase  of  the  data.  Also 
supporting  this  statement  is  the  general  agreement  of  form  drag  coefficients 
calculated  using  the  one  harmonic  approximation  to  those  calculated  using  the 
three  harmonic  fit  (which  follows  the  data  very  closely).  The  presence  of 
nonsymmetric  pressure  profiles  only  in  the  regions  where  separation  is 
predicted  by  theory  tends  to  support  these  predictions. 
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TABLE  A.l  PRESSURE  DATA  2ad/X  =  0.03125,  ct+=  0.0016 


Date:  June  3,1988 

Fluid:  Water  Temp:  25°  C  Density:  0.997  g/cm^  u:  .0.00896  cm^/s 
u*:  6.871  cm/s  a+:  60.8  Re:  41,600 
slope:  0.3906  intercept:  -0.3678 

Harmonic  1:  amplitudes  33.6  phases  +167° 

Harmonic  2:  amplitudes  0.97  phases  -177® 

Harmonic  3:  amplitudes  0.52  phases  +36® 


3.125E-02 

6.250E-02 

9.375E-02 

1.250E-01 

1.563E-01 

1.875E-OI 

2.188E-OI 

2.8I3E-01 

3.125E-01 

3.438E-01 

3.750E-01 

4.063E-01 

4.375E-01 

4.688E-0I 

5.000E-01 

5.313E-01 

5.625E-01 

5.938E-01 

6.250E-01 

6.563E-01 

6.875E-01 

7.188E-01 

7.500E-01 

7.813E-01 

8.125E-01 

8.438E-01 

8.750E-01 

9.063E-01 

9.375E-01 

9.688E-01 

1.000E+00 


-3.875E-01 

1.104E+00 

3.440E+00 

9.284E+00 

1.282E+01 

1.847E+01 

2.369E+01 

3.644E+01 

4.467E+01 

4.869E+01 

5.507E+0I 

5.858E+01 

6.218E+01 

6.541E+01 

6.902E+01 

6.958E+01 

6.763E+01 

6.726E+01 

6.456E+01 

6.204E+01 

5.775E+0I 

5.229E+01 

4.504E+01 

3.531E+01 

2.992E+01 

2.466E+01 

1.914E+01 

1.268E+01 

8.239E+00 

7.451E+00 

4.076E+00 


9.198E-01 

9.890E-01 

1.099E+00 

1.379E+00 

1.548E+00 

1.819E+00 

2.260E+00 

2.680E+00 

3.076E+00 

3.268E+00 

3.574E+00 

3.741E+00 

3.913E+00 

4.066E+00 

4.238E+00 

4.262E+00 

4.164E+00 

4.I44E+00 

4.010E+00 

3.884E+00 

3.673E+00 

3.406E+00 

3.05IE+00 

2.576E+00 

2.312E+00 

2.054E+00 

1.784E+00 

1.467E+00 

1.24SE+00 

1.208E+00 

I.04IE+00 


65 


TABLE  A.2  PRESSURE  DATA  2*dA  =  0.03125,  o+=  0.0021 
Date:  May  29,1988 

Fluid:  Water  Temp:  25°  C  Density:  0.997  g/cm^  »:  .0.00896  cmfys 
u*:  5.216  cm/s  a+:  46.2  Re:  30,400 
slope:  0.1670  intercept:  -0.0589 

Harmonic  1:  amplitudes  28.9  phases  +171° 

Harmonic  2:  amplitudes  0.82  phases  -132° 

Harmonic  3:  amplitudes  0.45  phases  -128° 

x  value/lambda  Differential  Pressure _ Average  Voltage 


3.125E-02 

5.374E+00 

3.775E-01 

6.250E-02 

6.635E+00 

4.557E-01 

9.375E-02 

8.522E+00 

5.748E-01 

1.250E-01 

1.251E+01 

8.312E-01 

1.563E-01 

1.672E+01 

1.102E+00 

1.875E-01 

2.142E+01 

1.405E+00 

2. 188E-01 

2.492E+01 

1.629E+00 

2.500E-01 

3.057E+01 

1.994E+00 

2.813E-01 

3.512E401 

2.288E+00 

3.125E-01 

4.048E+01 

2.634E+00 

3.438E-01 

4.657E+01 

3.302E+00 

4.063E-01 

5.642E+01 

3.663E+00 

4.375E-01 

6.004B+01 

3.895E+00 

4.688E-01 

6.068E+01 

3.933E+00 

5.000E-01 

6.537E+01 

4.235E+00 

5.313E-01 

6.216E+01 

4.021E+00 

5.625E-01 

6.230E+01 

4.026E+00 

5.938E-01 

6.231E+01 

4.022E400 

6.250E-01 

5.893E+01 

3.798E+00 

6.563E-01 

5.654E+01 

3.638E+00 

6.875E-01 

5.100E+01 

3.271E+00 

7.188E-01 

4.697E+01 

3.004E+00 

7.500E-01 

4.187E+01 

2666E+00 

7.813E-01 

3.493E+01 

2.209E+00 

8.125E-01 

2.900E+01 

1.817E+00 

8.438E-01 

2.342E+01 

1.448E+00 

8.750E-01 

1.802E+01 

1  092E+00 

9.063E-01 

1.521E+01 

9.039E-01 

9.375E-01 

1.026E+01 

5.763E-01 

9.688E-01 

8.156E+00 

4.349E-01 

1.000E+00 

6.312E+00 

3.102E-01 
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TABLE  A.3  PRESSURE  DATA  htfk  =  0.03125,  a+=  0.0029 


Date:  June  23,1988 

Fluid:  Water  Temp:  25°  C  Density:  0.997  g/cm^  t>:  .0.00896  cm^/s 
u*:  3.87  cm/s  a+:  34.3  Re:  21,600 

slope:  0.0922  intercept:  -0.0530 

Harmonic  1:  amplitudes  26.3  phases  +168° 

Harmonic  2:  amplitudes  1.47  phases  -122° 

Harmonic  3:  amplitudes  1.06  phases  -88° 

&  valuc/lambda _ Differential  Pressure _ Average.  Valine 


3.125E-02 

1.149E+00 

6.599E-01 

6.250E-02 

3.035E400 

7.812E-01 

9.375E-02 

4.341E-KX) 

8.639E-01 

1.250E-01 

7.512E+00 

1.071E+00 

I.563E-01 

1.190E+01 

1.359E+00 

1.875E-01 

1.432E+01 

1.516E+00 

2.188E-01 

1.729E+01 

1.709E+00 

2.500E-01 

2.273E+01 

2.067E+00 

2.813E-01 

2.497E+01 

2.212E+00 

3.125E-01 

3.142E+01 

2.638E+00 

3.438E-01 

3.557E+01 

2.909E+00 

3.750E-01 

4.265E+01 

3.377E+00 

4.063E-01 

4.724E+01 

3.678E+00 

4.375E-01 

4.958E+01 

3.830E+00 

4.688E-01 

5.076E+01 

3.904E+00 

5.000E-01 

5.240E+01 

4.009E+00 

5.313E-01 

5.140E+01 

3.938E+00 

5.625E-01 

5.276E+01 

4.024E+00 

5.938E-01 

5.139E+01 

3.929E+00 

6.250E-01 

5.187E+01 

3.957E+00 

6.563E-01 

4.869E+01 

3.741E+00 

6.875E-01 

4.431E+01 

3.445E+00 

7.188E-01 

4.012E+01 

3.J62E+00 

7.500E-01 

3.511E+01 

2.824E+00 

7.813E-OI 

2.757E+01 

2.317E+00 

8.125E-01 

2.251E+01 

1.977E+00 

8.438E-01 

1.709E+01 

1.612E+00 

8.750E-01 

1.230E+01 

I.289E+00 

9.063E-01 

7.020E+00 

9.325E-01 

9.375E-01 

3. 906 E +00 

7.209E-01 

9.688E-01 

1.841E+00 

5.792E-01 

1.000E+00 

8.668E-01 

5.101E-01 
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TABLE  A.4  PRESSURE  DATA  2ad fk  =  0.03125,  a+=  0.0040 


Date:  June  25,1988 

Fluid:  Water  Temp:  25°  C  Density:  0.997  g/cm^  o:  .0.00896  cm^/s 
a*:  2.80  cm/s  a+:  24.8  Re:  14,900 
slope:  0.0653  intercept:  -0. 1 1 24 

Harmonic  1:  amplitudes  20.3  phase=  4-169“ 

Harmonic  2:  amplitudes  2.23  phases  -109° 

Harmonic  3:  amplitudes  1.04  phases  -124° 

x  value/lambda  Differential  Pressure — Average  Voltage 


3.125E-02 

7.798E+00 

2.074E+00 

6.250E-02 

8.418E+00 

2.100E+00 

9.375E-02 

1.008E+01 

2.177E+00 

1.250E-01 

1.215E+01 

2.272E+00 

1.563E-01 

1.574E+01 

2.440E+00 

1.875E-01 

1.798E+01 

2. 5 44 E +00 

2.188E-01 

1.927E+0I 

2.602E+00 

2.500E-01 

2.368E+01 

2.808E+00 

2.813E-01 

2.485E+01 

2.862E+00 

3.125E-01 

2.796E+01 

3.006E+00 

3.438E-01 

3.249E+01 

3.219E+00 

3.750E-01 

3.526E+01 

3.348E+00 

4.063E-01 

3.974E+01 

3.558E+00 

4.375E-01 

4.330E+01 

3.725E+00 

4.688E-01 

4.489E+01 

3.797E+00 

5.000E-01 

4.733E+01 

3.910E+00 

5.313E-01 

4.638E+01 

3.862E+00 

5.625E-01 

4.478E+01 

3.783E+00 

5.938E-01 

4.565E+01 

3.822E+00 

6.250E-01 

4.518E+01 

3.796E+00 

6.563E-01 

4.357E+01 

3.716E+00 

6.875E-01 

4.049E+01 

3.567E+00 

7.I88E-01 

3.584E+01 

3.342E+00 

7.500E-01 

3.140E+01 

3.128E+00 

7.813E-01 

2.690E+0I 

2.911E+00 

8.125E-01 

2.121E401 

2.637E+00 

8.438E-01 

1.757E+01 

2.461E+00 

8.750E-01 

1.419E+01 

2.296E+00 

9.063E-01 

1.132E+01 

2.157E+00 

9.375E-01 

5.008E+00 

1.853E+00 

9.688E-0I 

4.900E+00 

1.845E+00 

1.000E+00 

4.158E+00 

1.807E+00 

TABLE  A.5  PRESSURE  DATA  2ad/X  =  0.03125,  o+=  0.011 


Date:  July  3,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.20  g/cm^  o:  .0.0805  cm' 
u*:  9.18  cm/s  a+:  9.06  Re:  4690 

slope:  0.2515  intercept:  -0.0653 

Harmonic  1:  amplitudes  12.7  phases  +150° 

Harmonic  2:  amplitudes  1.83  phases  -161° 

Harmonic  3:  amplitudes  0.46  phases  +52° 

x  value/lambda _ Differential  Pressure _ Average  Voltage 


3.125E-02 

4.488E-01 

3.221  E-01 

6.250E-02 

8.096E-01 

3.702E-01 

9.375E-02 

1.889E+00 

5.347E-01 

1.250E-01 

3.30IE+00 

7.528E-01 

1.563E-01 

5.082E+00 

1.031E+00 

1.875E-01 

6.851E+00 

1.307E+00 

2.188E-01 

8.649E+00 

1.588E+00 

2.500E-01 

1.062E+01 

1.896E+00 

2.813E-01 

1 .278E+01 

2.235E+00 

3.125E-01 

1.452E+01 

2.507E+00 

3.438E-01 

1.622E+01 

2.771E+00 

3.750E-01 

1.779E+01 

3.01 6E+00 

4.063E-01 

1.950E+01 

3.282E+00 

4.375E-01 

2. 102E+01 

3.518E+00 

4.688E-0 1 

2.247E+01 

3.742E+00 

5.000E-01 

2.417E+01 

4.007E+00 

5.313E-01 

2.542E+01 

4.200E+00 

5.625E-01 

2.699E+01 

4.442E+00 

5.938E-01 

2.791E+01 

4.582E+00 

6.250E-01 

2.839E+01 

4.649E+00 

6.563E-01 

2.824E+01 

4.615E+00 

6.875E-01 

2.715E+01 

4.428E+00 

7.188E-01 

2.526E+01 

4.1 12E+00 

7.500E-01 

2.303E+01 

3.741E+00 

7.813E-01 

1.994E+0I 

3.231E+00 

8.125E-01 

I.697E+01 

2.739E+00 

8.438E-01 

1.402E+01 

2.252E+00 

8.750E-01 

1.125E+01 

1.793E+00 

9.063E-01 

8.654E+00 

1.363E+00 

9.375E-01 

6.230E+00 

9.600E-01 

9.688E-01 

4.687E+00 

6.999E-01 

1.000E+00 

3.735E+00 

5.357E-01 

TABLE  A.6  PRESSURE  DATA  2ad fk  ■  0.03125,  a+=  0.015 


Dale:  July  1,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.20  g/cm^  v:  .0.0805  cm' 
u*:  6.661  cm/s  a+:  6.57  Re:  3250 

slope:  0.1402  intercept:  -0.0057 

Harmonic  1:  amplitudes  9.8  phases  +149° 

Harmonic  2:  amplitudes  1.8  phases  -161° 

Harmonic  3:  amplitudes  0.08  phases  -1° 


A 


value/lamhda  Differential  Pressure  Average  Voltage 


3.125E-02 

3.438E+00 

5.560E-01 

6.250E-02 

3.734E+00 

5.917E-01 

9.375E-02 

4.310E+00 

6.699E-01 

1.250E-01 

5.559E+00 

8.509E-01 

1.563E-01 

6.9S9E+00 

1.055E+00 

1.875E-01 

8.416E+00 

1.268E+00 

2.188E-01 

l.OOOE+OI 

1 .500E+00 

2.500E-01 

1.157E+01 

1.730E+00 

2.813E-01 

1.191E+01 

1.773E+00 

3.125E-01 

1.376E+01 

2.045E+00 

3.438E-01 

1.502E+01 

2.228E+00 

3.750E-01 

1.668E+0I 

2.472E+00 

4.063E-01 

1.767E+01 

2.614E+00 

4.375E-01 

1.903E+01 

2.812E+00 

4.688E-01 

2.023E+01 

2.985E+00 

5.000E-01 

2.170E+01 

3.200E+00 

5.313E-01 

2.165E+01 

3.183E+00 

5.625E-01 

2.278E+01 

3.346E+00 

5.938E-01 

2.429E-KU 

3.567E+00 

6.250E-01 

2.477E+01 

3.631E+00 

6.563E-01 

2.463E+01 

3.599E+00 

6.875E-01 

2.418E+01 

3.521E+00 

7.188E-01 

2.326E+01 

3.371E+00 

7.500E-01 

2.1 13E+01 

3.036E+00 

7.813E-01 

1.870E+01 

2.655E+00 

8.125E-01 

1.637E+01 

2.289E+00 

8.438E-01 

1.378E+01 

1.884E+00 

8.750E-01 

1.146E+01 

1.520E+00 

9.063E-01 

9.263E+00 

I.174E+00 

9.375E-01 

7.174E+00 

8.457E-01 

9.688E-01 

5.883E+00 

6.388E-01 

1.000E+00 

4.735E+00 

4.538E-01 
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TABLE  A.7  PRESSURE  DATA  2a<iA  =  0.03125,  a+=  0.015 
Date:  July  5,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.20  g/cm^  .0.0805  cm^/s 
u*:  6.60  cm/s  a+:  6.51  Re:  3220 


slope:  0.1032 

intercept:  -0.1766 

Harmonic 

1: 

amplitudes  10.2 

phases  +144° 

Harmonic 

2: 

amplitudes  1.38 

phases  +173° 

Harmonic 

3: 

amplitudes  0.53 

phases  +60° 

&__yalus/lanibda _ Diffasntial _ Pressure _ Average  Voltage 


3.125E-02 

-3.931E+00 

8.964E-01 

6.250E-02 

-3.587E+00 

9.537E-01 

9.375E-02 

-2.766E-KK) 

1.108E+00 

1.250E-01 

-1.718E+00 

1.309E+00 

1.563E-01 

-3.719E-01 

1.571E+00 

1.875E-01 

1.122E+00 

1.862E+00 

2.188E-01 

2.664E+00 

2. 164E+00 

2.500E-01 

4.294E+00 

2.484E+00 

2.813E-01 

5.958E+00 

2.810E+00 

3.125E-01 

7.482E+00 

3.108E+00 

3.438E-01 

8.977E+00 

3.400E+00 

3.750E-01 

1.026E+01 

3.649E+00 

4.063E-01 

1.168E+01 

3.925E+00 

4.375E-01 

1.288E+01 

4.157E+00 

4.688E-01 

1.393E+01 

4.358E+00 

5.000E-01 

1.532E+01 

4.628E+00 

5.313E-01 

1.647E+01 

4.851E+00 

5.625E-01 

1.741E+01 

5.029E+00 

5.938E-01 

1.855E+01 

5.250E+00 

6.250E-01 

1.903E+O1 

5.333E+00 

6.563E-01 

1.940E+01 

5.396E+00 

6.875E-01 

1.795E+01 

5.088E+00 

7.188E-01 

1.748E+01 

4.979E4-00 

7.500E-01 

1.589E+01 

4.641E+00 

7.813E-01 

1.355E+01 

4.152E+00 

8.125E-01 

1.153E+01 

3.726E+00 

8.438E-01 

9.284E+00 

3.256E+00 

8.750E-01 

7.083E+00 

2.794E+00 

9.063E-01 

4.947E+00 

2.345E+00 

9.375E-01 

3.1 13E+00 

1.953E+00 

9.688E-01 

1.787E+00 

1.674E+00 

1.000E+00 

1.035E+00 

1.508E+00 
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TABLE  A.8  PRESSURE  DATA  2a d/L  =  0.03125,  «+=  0.020 


Dale:  July  2,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1 .20  g/cm^  v>:  .0.0805  cm^/s 
u*:  5.02  cm/s  a+:  4.95  Re:  2360 

slope:  0.1384  intercept:  -0.2203 

Harmonic  1:  amplitudes  7.74  phases  +146° 

Harmonic  2:  amplitudes  1.57  phases  +174° 

Harmonic  3:  amplitudes  0.23  phases  +23° 

x value/lambda _ Differential  Pressure  Average  Voltage 


3.125E-02 

1.553E+00 

1.723E+00 

6.250E-02 

1.780E+00 

1 .737E+00 

9.375E-02 

2.502E+00 

1.795E+00 

1.250E-01 

3.619E+00 

1.888E+00 

1.563E-01 

5.079E+00 

2.01 1E+00 

1.875E-01 

5.998E+00 

2.086E+00 

2.188E-01 

7.486E+00 

2.212E+00 

2.500E-01 

8.823E+00 

2.324E+00 

2.813E-01 

9.890E+00 

2.412E+00 

3.125E-01 

1.061E+01 

2.470E+00 

3.438E-01 

1.156E+01 

2.548E+00 

3.750E-01 

1.254E+01 

2.629E+00 

4.063E-01 

1.372E+01 

2.727E+00 

4.375E-01 

1.462E+01 

2.801E+00 

4.688E-01 

1.551E+01 

2.873E+00 

5.000E-01 

1.646E+01 

2.952E+00 

5.313E-01 

1.744E+01 

3.032E+00 

5.625E-01 

1.715E+01 

3.001E+00 

5.938E-01 

1.834E+01 

3.100E+00 

6.250E-01 

1.850E+01 

3. 108E+00 

6.563E-01 

1.915E+01 

3.160E+00 

6.875E-01 

1.855E+01 

3.102E+00 

7.188E-01 

1.797E+01 

3.045E+00 

7.500E-01 

1.688E+01 

2.944E+00 

7.813E-01 

1.614E+01 

2.873E+00 

8.125E-01 

1.376E+01 

2.658E+00 

8.438E-01 

1.190E+01 

2.489E+00 

8.750E-01 

9.883E+00 

2.305E+00 

9.063E-01 

7.988E+00 

7.133E+00 

9.375E-01 

6. 1 19E+00 

1.963E+00 

9.688E-01 

4.770E+00 

1.838E+00 

1.000E+00 

4.061E+00 

1.770E+00 

I 


72 


i 


TABLE  A.  10  PRESSURE  DATA  2ad fK  =  0.0123,  a+=  0.00154 


Date:  September  3,1988 

Fluid:  Water  Temp:  25°  C  Density:  0.997  g/cm^  .0.00896  cm^/s 
u*:  7.205  cm/s  a+:  25.1  Re:  43.900 
slope:  0.1467  intercept:  -0.0244 

Harmonic  1:  amplitudes  12.59  phases  +172° 

Harmonic  2:  amplitudes  0.52  phases  +158° 

Harmonic  3:  amplitudes  0.22  phases  +85° 


3.125E-02 

2.347E+00 

4.904E-01 

6.250E-02 

3.548E+00 

6.519E-01 

9.375E-02 

4.462E+00 

7.725E-01 

1.250E-01 

5.770E+00 

9.491E-01 

1.563E-01 

8.167E+00 

1.280E+00 

1.875E-01 

9.09  IE +00 

1.402E+00 

2.188E-01 

1.234E+01 

1.855E+00 

2.500E-01 

1.441E+01 

2.139E+00 

2.813E-01 

1.772E+01 

2.599E+00 

3.125E-01 

1.994E+01 

2.905E+00 

3.438E-01 

2.218E+01 

3.215E+00 

3.750E-01 

2.383E+01 

3.440E+00 

4.063E-01 

2.489E+01 

3.582E+00 

4.375E-01 

2.646E+01 

3.795E+00 

4.688E-01 

2.674E+01 

3.825E+00 

5.000E-01 

2.842E+01 

4.055E+00 

5.3I3E-01 

2.870E+01 

4.086E+00 

5.625E-01 

2.814E+01 

3.998E+00 

5.938E-01 

2.665E+01 

3.777E+00 

6.250E-01 

2.659E+01 

3.759E+00 

6.563E-01 

2.487E+01 

3.507E+00 

6.875E-01 

2.288E+01 

3.215E+00 

7.188E-01 

2.022E+01 

2.828E+00 

7.500E-01 

1.798E+01 

2.501  E+00 

7.813E-01 

1.573E+01 

2.174E+00 

8.125E-01 

1.323E+01 

1.810E+00 

8.438E-01 

1.177E+01 

1.593E+00 

8.750E-01 

9.128E+00 

1.210E+00 

9.063E-01 

7.I30E+00 

9. 169E-01 

9.375E-01 

4.122E+00 

4.81  IE-01 

9.688E-01 

5.280E+00 

6.365E-01 

1.000E+00 

2.539E+00 

2.385E-01 

l 


I 
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TABLE  A.ll  PRESSURE  DATA  2a<iA  =  0.0123,  a+=  0.0020 

Dale:  August  31,1988 

Fluid:  Water  Temp:  25°  C  Density:  0.997  g/cm^  v:  .0.00896  cm 

u*:  5.483  cm/s  a+:  19.1  Re:  32,100 

slope:  0.0795  intercept:  -0.1349 

Harmonic  1:  amplitude=  10.85  phase=  +172° 

Harmonic  2:  amplitudes  0.47  phase=  -169'' 

Harmonic  3:  amplitudes  0.13  phases  +53° 

x  value/lambda  Differential  Pressure _ Average  Voltage 

3.125E-02 

2.605E-01 

1.550E+00 

6.250E-02 

1.035E+00 

1.656E+00 

9.375E-02 

1.785E+00 

1.758E+00 

1.250E-01 

2.727E+00 

1.889E+00 

1.563E-01 

4.383E+00 

2.126E+00 

1.875E-01 

5.167E400 

2.234E+00 

2.I88E-01 

7.596E+00 

2.586E+00 

2.500E-01 

9.600E+00 

2.875E+00 

2.813E-01 

1.292E+01 

3.360E+00 

3.125E-01 

1.505E+01 

3.668E+00 

3.438E-01 

1.649E+01 

3.873E+00 

3.750E-01 

1.78SE+01 

4.067E+00 

4.063E-01 

1.863E+01 

4.174E+00 

4.375E-01 

2.070E+01 

4.473E+00 

4.688E-01 

2.069E+01 

4.461  E+00 

5.000E-01 

2.232E+01 

4.695E+00 

5.313E-01 

2.167E+01 

4.590E+00 

5.625E-01 

2.100E+01 

4.480E+00 

5.938E-01 

1.975E+01 

4.283E+00 

6.250E-01 

2.007E+01 

4.323E+00 

6.563E-01 

1.878E+01 

4.121  E+00 

6.875E-01 

1.786E+01 

3.974E+00 

7.188E-01 

1.496E+01 

3.532E+00 

7.500E-01 

1.344E+01 

3.297E+00 

7.813E-01 

1.046E+01 

2.842E+00 

8.125E-01 

8.120E+00 

2.485E+00 

8.438E-01 

7.110E+00 

2.325E+00 

8.750E-01 

5.094E+00 

2.015E+00 

9.063E-01 

3.377E+00 

\  .749E+00 

9.375E-01 

2.061E+00 

1544E+00 

9.688E-01 

-1.331E-01 

1.207  E+00 

1.000E+00 

-1.191E+00 

1.040E+00 
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TABLE  A.  12  PRESSURE  DATA  2a<|A  =  0.0123,  a+-  0.0029 

Date:  September  2,1988 

Fluid:  Water  Temp:  25°  C  Density:  0.997  g/cm^  u:  .0.00896  cm 

u*:  3.843  cm/s  a+:  13.4  Re:  21,400 

slope:  0.0736  intercept:  -0.0319 

Harmonic  1:  amplitudes  8.98  phase=  +172° 

Harmonic  2:  amplitudes  0.74  phases  -175° 

Harmonic  3:  amplitudes  0.27  phases  -169° 

x  value/lambda  Differential  Pressure  Averase  Voltace 

3.125E-02 

3.793E+00 

7.335E-01 

6.250E-02 

5.064E+00 

8.307E-01 

9.375E-02 

6.166E+00 

9.143E-01 

1.250E-01 

6.705E+00 

9.525E-01 

1.563E-01 

7.861  E400 

1.040E+00 

1.875E-01 

8.996E+00 

1.I27E+00 

2.188E-01 

1.172E+01 

1.340E+00 

2.500E-01 

1.276E+01 

1.420E+00 

2.813E-01 

1.478E401 

1.577E+00 

3.125E-01 

1.631E+01 

1.695E+00 

3.438E-01 

1.782E401 

1.811E+00 

3.750E-0I 

1.868E+01 

1.875E+00 

4.063E-01 

1.865E+01 

1.868E+00 

4.375E-01 

2.120E401 

2.068E+00 

4.688E-0I 

2.166E+01 

2.100E+00 

5.000E-01 

2.250E+01 

2.162E+00 

5.313E-01 

2.319E401 

2.213E+00 

5.625E-01 

2.139E+01 

2.063E+00 

5.938E-01 

1.864E401 

1.836E+00 

6.563E-01 

2.008E+01 

1.942E+00 

6.875E-01 

1.930E4O1 

1.874E-KX) 

7.188E-01 

1.682E+01 

1 .669E+00 

7.500E-01 

1.539E+01 

1.549E+00 

7.8I3E-01 

1.356E+01 

1.397E+00 

8.125E-01 

I.144E+01 

1.221E+00 

8.438E-01 

1.072E+01 

1.158E+00 

8.750E-01 

8.800E+00 

9.986E-01 

9.063E-01 

6.463E+00 

8.053E-01 

9.375E-01 

5.283E+O0 

7.053E-01 

9.688E-01 

5.176E+00 

6.916E-01 

1.000E+00 

1.353E+00 

3.789E-01 

I 

I 
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TABLE  A.13  PRESSURE  DATA  2ad/X  =  0.0123,  a+=  0.0039 


Date:  September  14,  1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.122  g/cm^  m:  .0.0249  cm^/s 
u*:  7.84  cm/s  a+:  9.85  Re:  15,000 

slope:  0.2031  intercept:  -0.0584 


Harmonic  1 
Harmonic  2 
Harmonic  3 


amplitudes  9.78 
amplitudes  0.60 
amplitudes  0.10 


phases  +163° 
phases  +158° 

phases  +28° 


3.125E-02 

1.181E+00 

4.400E-01 

6.250E-02 

1 .930E+00 

5.335E-01 

9.375E-02 

2.635E400 

6.21  IE-01 

1.250E-01 

3.602E+00 

7.442E-01 

1.563E-01 

5.184E+00 

9.514E-01 

1.875E-01 

6.546E+00 

1.129E+00 

2.188E-01 

8.579E+00 

1.397E+00 

2.500E-01 

1.047E+01 

1.646E+00 

2.813E-01 

1.181E+01 

1.821E+00 

3.125E-01 

1.379E+01 

2.082E+00 

3.438E-01 

1.551E+01 

2.308E+00 

3.750E-01 

1.723E+01 

2.534E+00 

4.063E-01 

1.882E+01 

2.743E+00 

4.375E-G1 

2.041E+01 

2.950E+00 

4.688E-01 

2.117E+01 

3.045E+00 

5.000E-01 

2.254E-KJ1 

3.223E+00 

5.313E-01 

2.067E401 

2.960E+00 

5.625E-01 

2.127E+01 

3.033E+00 

5.938E-01 

2.094E+01 

2.979E+00 

6.250E-01 

2.110E+0I 

2.993E+00 

6.563E-01 

2.048E+01 

2.900E+00 

6.875E-01 

1.973E+01 

2.789E+00 

7.188E-01 

1.736E+01 

2.457E+00 

7.500E-01 

1.636E+01 

2.31 1E+00 

7.813E-01 

1.318E+01 

1.869E+00 

8.125E-01 

1.126E+01 

1.598E+00 

8.438E-01 

1.022E+01 

1.447E+00 

8.750E-01 

8.603E400 

1.219E+00 

9.063E-01 

6.634E+00 

9.415E-01 

9.375E-01 

4.679E+00 

6.660E-01 

9.688E-01 

3.568E+O0 

5.059E-01 

1.000E+00 

3.550E+00 

4.948E-01 
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TABLE  A.  14  PRESSURE  DATA  2adA  =  0.0123,  a+=  0.0059 

Date:  September  11,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.159  g/cm^  v:  .0.0406  cm^/j 
u*:  8.542  cm/s  a+:  6.57  Re:  9450 

slope:  0.1316  intercept:  -0.1415 

Harmonic  1:  amplitudes  6.79  phases  +162° 

Harmonic  2:  amplitudes  0.36  phases  -133° 

Harmonic  3:  amplitudes  0.08  phases  +12° 

x  value/lambda _ Differential  Pressure  Average  Voltage 

3.125E-02 

-2.197E-01 

8.402E-01 

6.250E-02 

-7.968E-02 

8.599E-01 

9.375E-02 

3.081E-01 

9.438E-01 

1.250E-01 

9.055E-01 

1.082E+00 

1.563E-01 

1 .883E+00 

1.318E+00 

1.875E-01 

2.526E+00 

1.468E+00 

2.188E-01 

3.327E-KX) 

1.658E+00 

2.500E-01 

4.697E+00 

1  996E+00 

2.813E-01 

6.167E+00 

2.360E+00 

3.125E-01 

7.467E+00 

2.679E+00 

3.438E-01 

8.599E+00 

2.955E+00 

3.750E-0I 

9.632E+00 

3.206E+00 

4.063E-01 

1.020E+01 

3.337E+00 

4.375E-01 

1.161E+01 

3.686E+00 

4.688E-01 

1.208E+01 

3.790E+00 

5.000E-01 

1.300E+01 

4.010E+00 

5.313E-01 

1.319E401 

4.044E+00 

5.625E-01 

1.344E+01 

4.091E+00 

5.938E-01 

1.280E+01 

3.911E+00 

6.250E-01 

1.316E+01 

3.986E+00 

6.563E-01 

1.261E+01 

3.827E+00 

6.875E-01 

1.166E+01 

3.567E+00 

7.188E-01 

1.026E401 

3.188E+00 

7.500E-0! 

7.562E+00 

2.458E+00 

8.125E-01 

6.1 19E+00 

2.068E+00 

8.438E-01 

5.058E-KX) 

1.778E+00 

8.750E-01 

3.380E+00 

1.328E+00 

9.063E-01 

2.169E+00 

9.983E-01 

9.375E-01 

9.476E-01 

6.663E-01 

9.688E-01 

5.270E-01 

5.412E-01 

1.000E+00 

3.347E-01 

4.750E-01 

I 


78 


TABLE  A.15  PRESSURE  DATA  2a<jA  =  0.0123,  o+=  0.0079 

Date:  September  10,1988 

Fluid:  Dextrpse  Temp:  25°  C  Density:  1.159  g/cm^  .0.0406  cm 

u*:  6.35  cm/s  a+:  4.88  Re:  6730 

slope:  0.0959  intercept:  -0.0573 

Harmonic  1:  amplitudes  5.72  phase=  +157° 

Harmonic  2:  amplitudes  0.52  phases  -127° 

Harmonic  3:  amplitudes  0.32  phases  +44° 

x  value/lambda _ Differential _ Pressure _ Average  Voltage 

3.125E-02 

-7.393E-01 

4.40  IE-01 

6.250E-02 

-5.819E-01 

4.585E-01 

9.375E-02 

-1.911E-01 

5.225E-01 

1.250E-01 

-1.650E-01 

5. 152E-01 

1.563E-01 

7.121E-01 

6.744E-01 

1.875E-01 

1.707E+00 

8.565E-01 

2.188E-01 

2.297E+00 

9.596E-01 

2.500E-01 

2.973E+00 

1.080E+00 

2.813E-01 

4.409E+00 

1.348E+00 

3.125E-01 

5.515E+00 

1.552E+00 

3.438E-01 

6.483E+00 

1.729E+00 

3.750E-01 

6.972E+00 

1.8I2E+00 

4.063E-01 

7.627E+00 

1.928E+00 

4.375E-01 

8.600E+00 

2.I06E+00 

4.688E-01 

8.572E+00 

2.088E+00 

5.000E-01 

9.5I8E+00 

2.260E+00 

5.313E-01 

1.087E+01 

2.512E+00 

5.625E-01 

1.127E+01 

2.579E+00 

5.938E-01 

1.060E+01 

2.434E+00 

6.250E-01 

1.121E+01 

2.543E+00 

6.563E-01 

1.067E+01 

2.424E+00 

6.875E-01 

1.013E+01 

2.305E+00 

7.188E-01 

8.386E+00 

I.952E+00 

7.500E-01 

6.864E+00 

1.642E+00 

7.813E-01 

6.120E+00 

1 .484E+00 

8.125E-01 

5.104E+00 

1.273E+00 

8.438E-01 

4.014E+00 

I.047E+00 

8.750E-01 

2.491E+00 

7.369E-01 

9.063E-01 

1.356E+00 

5.024E-01 

9.375E-01 

5.421E+00 

1.285E+00 

9.688E-01 

7.481E-02 

2.270E-01 

1.000E+00 

-9.929E-02 

1.805E-01 

TABLE  A.  16  PRESSURE  DATA  2a<i/X  =  0.0123,  «+=  0.0099 


Date:  September  8,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.205  g/cm^  u:  .0.0740  cafi/s 
u*:  9.23  cm/s  a+:  3.90  Re:  5200 

slope:  0.1339  intercept:  -0.1560 

Harmonic  1:  amplitudes  4.69  phase=  +162° 

Harmonic  2:  amplitudes  0.38  phase=  +175° 

Harmonic  3:  amplitudes  0.17  phases  +10° 

x  value/lambda  Differential  Pressure  Average  Voltage 

3.125E-02 

1.426E+00 

1.585E+00 

6.250E-02 

1 .682E+00 

1.644E+00 

9.375E-02 

1.900E+00 

1.692E+00 

1.250E-01 

2.161E+00 

1.753E+00 

1.563E-01 

2.538E+00 

1.849E+00 

1.875E-01 

3.265E+00 

2.054E+00 

2.188E-01 

3.967E+00 

2.251E400 

2.500E-01 

5.300E+00 

2.642E+00 

2.813E-01 

6.465 E-KX) 

2.981E+00 

3.125E-01 

7.466E+00 

3.270E+00 

3.438E-01 

7.762E+00 

3.342E+00 

3.750E-01 

8.794E+00 

3  640E+00 

4.063E-01 

9.257E+00 

3.764E+00 

4.375E-01 

9.497E+00 

3.8I8E+00 

4.688E-01 

9.684E+00 

3.856E+00 

5.000E-01 

1.012E+01 

3.970E+00 

5.313E-01 

1.061E+01 

4. 103E+00 

5.625E-0I 

1.103E+0I 

4.212E+00 

5.938E-01 

9.987E+00 

3.871E+00 

6.250E-01 

1.099E+01 

4.161E+00 

6.563E-01 

1.027E+01 

3.919E+00 

6.875E-01 

9.786E+00 

3.751E+00 

7.188E-01 

8.874E+00 

3.450E+00 

7.500E-01 

7.978E+00 

3.154E+00 

7.813E-01 

7.213E+00 

2.899E+00 

8.125E-01 

6.539E-KX) 

2.672E+00 

8.438E-01 

5.367E-KX) 

2.291  E+00 

8.750E-01 

4.367E+00 

1.963E+00 

9.063E-01 

3.412E+00 

1 .649 E+00 

9.375E-01 

2.293E+00 

1.285E+00 

9.688E-01 

2.012E+00 

1.179E+00 

1.000E+00 

1.71 1E+00 

1.067E+00 
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TABLE  A.  17  PRESSURE  DATA  2ad/X  =  0.0123.  o+=  0.0101 


Dale:  September  7,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.205  g/cm^  .0.0740  cm 

u*:  9.102  cm/s  a+:  3.84  Re:  5120 

slope:  0.1109  intercept:  -0.0676 

Harmonic  1:  amplitudes  4.45  phases  +156° 

Harmonic  2:  amplitudes  0.28  phases  -147° 

Harmonic  3:  amplitudes  6.13  phases  +121° 

x  value/lambda  Differential _ Pressure  Average _ Voltage 

3.125E-02 

5.658E-01 

7.914E-01 

6.250E-02 

6.825E-01 

8.107E-01 

9.375E-02 

1.056E+00 

9.229E-01 

1.250E-01 

1.422E+00 

1.033E+00 

1.563E-01 

2.185E400 

1.286E+00 

1.875E-01 

2.844E+00 

1.501E+00 

2.188E-01 

3.563E+00 

1.739E+00 

2.500E-01 

4.271E-KX) 

1.972E+00 

2.813E-01 

4.388E+00 

1.991E+00 

3.125E-01 

5.278E+00 

2.291E+00 

3.438E-01 

6.191E+00 

2.598E+00 

3.750E-01 

6.890E+00 

2.828E+00 

4.063E-01 

7.640E+00 

3.077E+00 

4.375E-01 

8.410E+00 

3.333E+00 

4.688E-01 

8.972E+00 

3.514E+00 

5.000E-01 

9.479E400 

3.674E+00 

5.313E-01 

9.746E+00 

3.748E+00 

5.625E-01 

9.744E+00 

3.724E+00 

5.938E-01 

8.668E+00 

3.312E+00 

6.250E-01 

9.849E+00 

3.716E+00 

6.563E-01 

9.675E+00 

3.630E+00 

6.875E-01 

9.100E+00 

3.399E+00 

7.188E-01 

8.264E+00 

3.073E+00 

7.500E-01 

7.231E+00 

2.676E+00 

7.813E-01 

6.605E+00 

2.427E+00 

8.I25E-01 

5.748E+00 

2.094E+00 

8.438E-01 

4.857E-KX) 

1.748E+00 

8.750E-01 

3.754E+00 

1.326E+00 

9.063E-01 

2.987E+00 

1 .025E+00 

9.375E-01 

2.253E+00 

7.359E-01 

9.688E-01 

1 .805E+00 

5.510E-01 

1.000E+00 

1.621E+00 

4.614E-01 
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TABLE  A.  18  PRESSURE  DATA  2*dA  =  0.0123,  a+=  0.00148 

Date:  September  6,1988 

Fluid:  Dextrose  Temp:  25°  C  Density:  1.205  g/cm^  o:  .0.0740  cm^/ 
u*:  6.202  cm/s  a+:  2.62  Re:  3300 

slope:  0.0762  intercept:  -0.1094 

Harmonic  1:  amplitudes  3.72  phases  +164° 

Harmonic  2:  amplitudes  0.45  phases  -179° 

Harmonic  3:  amplitudes  0.11  phases  -134° 

x  value/lambda  Differential  Pressure _ Average _ Voltage 

3.125E-02 

3.159E-01 

1.497E+00 

6.250E-02 

8.044E-01 

1.601E+00 

9.375E-02 

9.874E-0I 

1.63 1  E+00 

1.250E-01 

9.791E-01 

1.613E+00 

1.563E-01 

1.581E+00 

1.745E+00 

1.875E-01 

2.294E+00 

1.904E+00 

2.188E-01 

2.981  E+00 

2.056E+00 

2.500E-01 

3.887E+00 

2.262E+00 

2.813E-01 

4.585E+00 

2.418E+00 

3.125E-01 

5.074E+00 

2.522E+00 

3.438E-01 

5.431E+00 

2.594E+00 

3.750E-01 

6.501E+00 

2.824E+00 

4.375E-01 

7.1 15E+00 

2.959E+00 

4.688E-01 

7.292E400 

2.987E+00 

5.000E-01 

7.9I8E+00 

3.125E+00 

5.313E-01 

7.552E+00 

3.020E+00 

5.625E-01 

7.638E+00 

3.025E+00 

5.938E-01 

4.741E+00 

2.301  E+00 

6.250E-01 

7.476E+00 

2.954E+00 

6.563E-01 

7.030E+00 

2.829E+00 

6.875E-01 

7.414E+00 

2.908E+00 

7.188E-01 

6.462E+00 

2.659E+00 

7.500E-01 

4.206E+00 

2.092E+00 

7.813E-01 

4.960E+00 

2.261  E+00 

8.125E-01 

4.225E+00 

2.066E+00 

8.438E-01 

3.824E+00 

1.952E+00 

8.750E-01 

2.597E+00 

1.636E+00 

9.063E-01 

1.903E+00 

1. 451 E+00 

9.375E-01 

7.825E-01 

1.1 61  E+00 

9.688E-01 

4.447E-01 

1.063E+00 

1.000E+00 

-1.913E-02 

9.342E-01 
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APPENDIX  B 

Dala  Analysis  Programs 


I 

I 
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C*********************************************************** 
£**********************************************************♦ 
c****«,*************PROORAM  PRESS.FOR*******  *******  ********** 

C*  This  program  reads  the  binary  file  created  by  the  Metra- 

C*  Byte  Analog  to  Digital  conversion  board  using  program 

C*  Stream  16.  There  are  two  files  necessary  to  run  this  pro 
C*  gram.  The  first  file,  for008.dat,  should  contain:  the 
C*  number  of  points  per  file,  value  of  intercept  and  slope 
C*  from  pressure  transducer  calibration,  density  of  fluid, 

C*  friction  velocity  and  a  listing  of  the  following  (at 
C*  each  x  value)  :  x  value(inches),  data  filename,  raw  data 
C*  output  index.  The  second  file,  for009.dat,  should  con- 
C*  tain:  raw  output  filename  and  the  number  of  points  in 
C*  that  file(at  each  x  value,  if  necessary),  and  the  out- 
C*  put  file  name  for  the  plotting  file. 

C*  The  program  will  convert  the  binary  data  to  ascii 

C*  form  with  units  of  voltage.  The  dimensionless  pressure  is  calculated 

C*  and  plotting  files  are  then  constructed  . 

C*  By  Dave  Courtemanche  3/2/88 

£********»*•************************************************ 

£*********************************************************** 


I 


I 

I 

I 


character*12  filenamel(32),  routname 
character*  12  profilename,  fitname 
integer*4  i,  j,  k,  ii,  jj,  kk,  ij,  nleft,  nrec 
integer*4  ipts,  nnleft,  nnrec 
integer*4  nopts,  nnopts,  numpts,  ndiff 
integer*2  buffer(256) 
real  xval(32),  value(30000),  val 
real  row,  ustar,  norm 

real  delp(lOOOO),  totp.totv,  a,  b,  avgp(32),avgv(32) 
real  int(32),  skew(32),  flat(32),  delpp(32),  intp(32) 

c*********************************************************** 
C*  nnopts  is  the  #  of  points  per  file  * 

£*****•*•**************************•*******•**********»***** 

read(8,*)  nnopts,  a,  b,  row,  ustar 
nnleft  =  mod(nnopts,256) 
nnrec=  nnopts/256 
norm  =  row*ustar*ustar/2486.2 

do  10  ij=  1 ,  32 

£*********************************************************** 

C*  xval(ij)  is  the  x  value(inches).  (This  will  be  con-* 

C*  verted  to  x/lambda.)  Filename  1  (ij)  is  the  raw  press* 

C*  data  filename  the  this  x  location.  Irawout  (ij)  * 

C*  should  be  1  if  you  want  to  output  raw  data  for  this  x 

£******»»•*****************************»******************** 
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read(8,3000)  xval(ij),  filenamel(ij) 
c*  @,  Irawout 

3000  format  (f8.6,  ail,  4i3) 
xval(ij)  =  xval(ij)  /  2.00 

open(unit=10,file=filcnamei(ij),form=,binary', 
<§>slatus=’oId’,recl=I6,  access='direct') 


c***************************************»*****< 
C*  This  part  of  program  unpacks  data  as  voltage,  converts 
C*  it  to  pressure(inches  of  water),  and  calculates  the 
C*  average.  * 

C*  ******* ***************************** •******•< 

kj=0 

nopts=nnopts 
nleft=nnleft 
nrec=nnrec  +1 
kk=0 
totp=0 
totv=0 
int(ij)=0 
skcw(ij)=0 
flat(ij)=0 

do  20  k=l,  nrec 
if(k  .It.  nrec)  then 
read  (10)  buffer 

do  30  ii  =  1.256 

call  unpack  (buffer(ii),  ichan,  val  ) 
kk  =  kk+  1 
value(kk)  =  val 
de!p(kk)  =  a  +  b  *  value(kk) 
totp  =  totp  +  delp(kk) 
totv  =  totv  +  value(kk) 

30  continue 

else 

read  (10)  buffer 
do  40  ii  =  1  jnleft 

call  unpack(buffer(ii),  ichan,  val) 
kk  =  kk  +  1 
value(kk)  =  val 
delp(kk)  =  a  +  b  *  value(kk) 
totp  =  totp  +  delp(kk) 
totv  =  totv  +  vaiue(kk) 


I 


40  continue 
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endif 

20  continue 

rewind  (10) 

av8P('j)  =  totp/nopts 
avgv(ij)  =  totv/nopts 

C*  Outputting  raw  voltage  data  (if  desired), 

c*  if  (irawout  .ne.l  )  goto  3 
c*  read(9,3400)  routname,  irawnum 

c*3400  format  (al2,i4) 

c*  open(unit=27,  file=routname,form='formatted',status= 

c*  (a) 'new') 

c*  do  50  i  =  1,  irawnum 

c*  write(27,3500)  value(i) 

c*  50  continue 

c*3500  format(lpel  1.3) 

c*  3  continue 

C*  CALCULATING  THE  HIGHER  MOMENTS************************* 
do  15  kk=l,nopts 

int(ij)»int(ij)+(delp(kk)*avgp(ij))**2 

skew(ij)=skew(ij)+(delp(kk)-avgp(ij))**3 

flat(ij)=flat(ij)+(delp(kk)-avgp(ij))**4 

15  continue 

int(ij)  =  sqrt(int(ij)/nopts) 

skew(ij)  =  skew(ij)/(nopts  *  int(ij)**3) 

flat(ij)  =  flat(ij)/(nopts  *  int(ij)**4) 

10  continue 

C*  Subtracting  mean  pressure  gradient 

do  37  kk=l,32 

avgp(kk)=  avgp(kk)  +  2.0305*norm*float(kk)/32. 

37  continue 

C*  Non-dimensionalizing  with  wall  parameters 

do  16  1=1,32 

delpp(l)=avgp(l)/norm 

intp(l)=int(l)/norm 

16  continue 

C***  outputting  plotting  files*********** 

read(9,3550)  profilename 
read(9,3550)  filname 
3550  format(al2) 

open(unit=50,file=profilename,form=’forrnatted'. 


@status='new') 

open(unit=5 1  ,file=fitname,form=formatted\ 
@status='new') 
write  (50,6350) 

6350  formal('xval/lambda  \2x,'  delp  \3x,’  deipp' 
@,3x,’  intp  \3xf'  avgv  ') 


do  100  11=1,32 

write(5 0,6400)  xval(ll),avgp(U),delpp(M), 
@intp(U),avgv(U) 

6400  format(6(lpel0.3,2x)) 

wrile(5 1,6500)  xval(ll),delpp(ll) 

6500  format(2(lpel0.3,2x)) 

100  coatinue 

17  stop 
end 

C*  Subroutine  used  to  unpack  data  from  binary  to  ascii 
subroutine  unpack  (ibin, ichan,  val) 
integer*2  ibin 
integer*  4  ichan,  ival 
real  val,  aval 

ichan  =  iand(ibin,  15) 
if  (ibin  .It.  2047)  then 
ival  =  iand(ibin,  65535) 

else 

ival  =  ibin 
endif 

ival  =  ival  /  16 
aval  =  ival 

val  =  (aval/4095.)  *5. 


return 

end 
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fit3u.for 

C*  This  program  reads  in  a  two  dimensional  array 
C*  It  then  fits  the  second  variable  of  the  array  to  the  first 
C*.. using  a  three  harmonic  Fourier  series. 

C*  A  plotting  file  is  output  which  shows  both  the  first  harmonic  fit 
C*  and  the  fit  using  all  three  harmonics, 
program  datfit 
dimension  a(20,20),b(20,l) 
di  mens  ion  amag(  1 0)  ,z(  1 00) ,  x(  1 00)  ,angd(  1 0)  ,ang(  1 0) 
dimension  work(20000)  ,rat ( 1 0)  ,xl(600)  ,f(600) 
dimension  ff(600) 

character*  1 2,  filename  1  ,filename2,filename3 
real  mean, pi 

print*,'  Enter  the  input  file  name' 
read(5 ,22)filename  1 

print*,'  Enter  the  number  of  points  (0  to  360  degrees)’ 
read*,np 

c  print*,'  three  harmonics  will  be  calculated  ' 

print*,'  Enter  harmonic  output  file  name' 
read(5,22)fiiename3 
22  format(al2) 

open(unit=40,file=filenamel  .form=’formatted',status=’old') 
open(unit=72,file=filename3,form='formatted',status=’new') 

nh  arm  o=3 
n=nharmo*2 
mo=n+l 

pi=4.0*atan(1.0) 

do  1  i=l  ,mo 
do  2  j=  1  ,mo 

2  a(ij)=0. 

1  continue 

do  3  k=l,mo 

3  b(k,l)=0. 

do  4  i=l,np 

read(40,433)x(i),z(i) 
x(i)=x(i)*2*pi 
4  continue 

do  5  m=l,np 
ii=0 

b(mo,  1  )=b(mo,  1  )+z(m) 
do  6  jj=l,n,2 
ii=ii+l 

b(jj,  1  )=b(jj,l  )+z(m)*sin(x(m)*ii) 
b(jj+1.0=b(jj+l,l)+z(m)*cos(x(m)*ii) 


mm=0 

do  7  kk=l  ,n,2 
mm=mm+l 

a(jj,kk)=a(jj,kk)+sin(x(m)*ii)*sin(x(m)*mni) 
a(jj,kk+ 1  )=a(jj^k+ 1  )+sinf  xf  m)*ii)*cosfxf  m)*mm) 
a(jj+l  ,kk)=a(jj+ 1  ,kk)+cos(x(m)*ii)*sinfxfm)*mm) 
a(jj+ 1  ,kk+ 1  )=a(jj+ 1  ,kk+ 1  )+cos(xf  m)*ii)*cosf  x(m)*mm) 

7  continue 

a(jj,mo)=a(jj,mo)+sin(x(tn)*ii) 
a(jj+ 1  ,mo)=a(jj+ 1 ,mo)+cos(xfm)*ii) 
a(mojj)=a(jjano) 
a(mo  jj+l)=a(jj+l  ,mo) 

6  continue 
5  continue 

a(mo,mo)=float(np) 
call  leqt  1  f(a,  1  ,mo,20,b,0, work.ier) 
print*,'  ier=  ',ier 

if=0 

do  8  ie=l,n,2 
if=if+l 

amag(if)=sqrt(b(ie,  1  )**2+b(ie+l ,  1  )**2) 
bb=-b(ie,l) 
bbb=b(ie+l,l) 

ang(if)=atan2(bb,bbb) 
angd(if)=ang(if)*360./(2.*pi) 
rat(if)=amag(if)/amag(  1 ) 

8  continue 

zavg=b(mo,l) 
print  lOO.zavg 

print  101,(ih,amag(ih),angd(ih),rat(ih),ih=l,3) 
write(72, 1 00)zavg 

write(72, 101  )(ih,amag(ih)  ,angd(ih),rat(ih),ih=  1 ,3) 

100  format(//,’  mean=',2x,fl2.5) 

101  format(  lx,'harmonic=  ',i2,3x,'amp=  ’,f9.5,3x,'phase=  ’,fl0.3, 
@  3x,'  amp.  ratio=  ’,f7.5) 

C**  This  section  generates  a  curve  of  the  fourier  fit 

print*,'  Include  the  mean?  l=yes  ' 
read*, reply 

mean=0.0 

ifereply.eq.l)mean=zavg 


do  800  i=l  ,250 

xKi)=noateei-l))/249.0 

sum=0.0 
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do  850  j=l,l 

sum=sum+amag(j)*cos(2.0*pi*j','xi(i)+ang(j)) 
850  continue 

f(i)=mean+sum 

sum  1=0.0 
do  860  j=l,3 

suml=suml+amag(j)*cos(2.0*pi*j*xl(i)+ang(j» 
860  continue 

ff(i)=mean+sum  1 
800  continue 

do  i=  1,250 

write(72,433)xl(i),f(i),ff(i) 
end  do 


433  format(3(lpe!0.3,2x)) 
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APPENDIX  C 
Sample  Calibration 


1 
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NOMENCLATURE 

A  van  Driest  parameter 

A  Average  van  Driest  parameter,  see  equation  (3.18) 

ad  amplitude  of  wave  surface  (cm) 

aj  amplitude  of  i*b  harmonic  of  non-dimensionalized  wave  induced 

surface  pressure 

a*  non-dimensional  wave  surface  amplitude 

bj  coefficients  in  orthogonal  transformation 

Cf  friction  factor 

h  average  channel  half  height  (cm) 

Ki,K2  non-dimensional  empirical  constants,  see  equation  (3.18) 

Kl  non-dimensional  empirical  lag  constant,  see  equation  (3.19) 

pd  surface  pressure  (inches  H2O) 

p  non-dimensional  surface  pressure  (p=pd /pu*2) 

p+  non-dimensional  pressure  gradient,  see  equation  (3.18) 

p+eff  effective  pressure  gradient,  see  equation  (3.19) 

Re  Reynolds  number,  based  on  bulk  velocity  and  channel  half 

height 

Rji  time-averaged  stress  tensor  (g/cm  s^) 

Sii  time-averaged  rate  of  strain  tensor  (1/s) 

u  streamwise  velocity  (cm/s) 

ub  bulk  velocity  (cm/s) 

u '  fluctuating  component  of  streamwise  velocity  (cm/s) 

u  *  friction  velocity  (cm/s) 

normal  velocity  (cm/s) 


v 


fluctuating  component  of  normal  wise  velocity  (cm/s) 
Cartesian  streamwise  coordinate  (cm) 


v’ 
x 

x+  non-dimensional  Cartesian  streamwise  coordinate 

y  Car(esian  normal  coordinate  (cm) 

y'  distance  to  solid  surface  (cm) 

yt  height  of  top  wall  of  channel  in  Cartesian  coordinates  (cm) 

Greek  letters 

a  wave  number  (1/cm) 

a*  dimensionless  wave  number 

q  curvilinear  normal  coordinate  (cm) 

tit  height  of  top  wall  in  curvilinear  coordinates  (cm) 

K  von  Karman  constant 

X  wavelength  of  wave  surface  (cm) 

0j  phase  angle  of  i1^  harmonic  of  non-dimensionalized  wave 

induced  surface  pressure 

v  kinematic  viscosity  (cm^/s) 

vt  turbulent  eddy  viscosity,  see  equations  (3.14-16)  (cm2/s) 

e  curvilinear  streamwise  coordinate  (cm) 

p  density  (g/cm^) 

x  local  shear  stress  (g/s^) 

V  stream  function,  see  equations  (3.9-10)  (cm^/s) 

©  vorticity,  see  equation  (3.8)  (1/s) 
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